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Abstract 

We investigate causality constraints on the time evolution of entanglement entropy after a global 
quench in relativistic theories. We first provide a general proof that the so-called tsunami velocity 
is bounded by the speed of light. We then generalize the free particle streaming model of [1] to 
general dimensions and to an arbitrary entanglement pattern of the initial state. In more than two 
spacetime dimensions the spread of entanglement in these models is highly sensitive to the initial 
entanglement pattern, but we are able to prove an upper bound on the normalized rate of growth 
of entanglement entropy, and hence the tsunami velocity. The bound is smaller than what one 
gets for quenches in holographic theories, which highlights the importance of interactions in the 
spread of entanglement in many-body systems. We propose an interacting model which we believe 
provides an upper bound on the spread of entanglement for interacting relativistic theories. In two 
spacetime dimensions with multiple intervals, this model and its variations are able to reproduce 
intricate results exhibited by holographic theories for a significant part of the parameter space. 
For higher dimensions, the model bounds the tsunami velocity at the speed of light. Finally, we 
construct a geometric model for entanglement propagation based on a tensor network construction 
for global quenches. 
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I. INTRODUCTION AND SUMMARY 

Understanding the evolution of quantum entanglement in non-equilibrium processes such 
as thermalization is a question of much interest. Entanglement could reveal quantum cor¬ 
relations not easily accessible by other observables such as thermodynamic quantities or 
correlation functions, and thermalization provides a dynamical setting to study the genera¬ 
tion and spread of entanglement between subsystems. 

A simplest physical context to probe this question is the evolution of entanglement en¬ 
tropy <S's(t) after a global quench in a conformal field theory (GFT), where one injects a 
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uniform energy density in a very short time interval at t = 0 and then lets the system evolve. 
Here E denotes the entangling surface, whose characteristic size R will be taken to be much 
larger than the inverse equilibrium temperature, i.e. R ^ 1/T. Recent studies of evolution 
of Sj:{t) for quench processes in (1 + 1) dimensions as well as in holographic systems of 
general dimensions have revealed a “universal” linear regime [1-3]^ 

ASj:{t) = VESeqAj:t, R>t>4q- (1-1) 

Here ASj:{t) is the change of <S's(t) from its value at t = 0, is the area of the entangling 
surface S, and Seq is the equilibrium entropy density. 4q ~ ^/T is the local equilibra¬ 
tion time. More explicitly, we expect for a region A whose size Ra is comparable to, but 
larger than i^q, that its entanglement entropy saturates at the thermal value after the local 
equilibration time, i.e. 

Sa = Seqi4, for Ra>^ and t» 4q (1.2) 

with Va the volume of A. 

In (1.1) ve is a constant of dimension velocity and depends on macroscopic properties of 
the state. With the speed of light set to c = 1, for (1 + 1) dimensions it was found [1] that 

VE = l, d = 2, (1.3) 


while for higher dimensional holographic systems at zero chemical potential [2, 3], 
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Equation (1.1) suggests a simple heuristic picture for the growth of entanglement: an 
entanglement wave propagates inward from the boundary of the entangled region, with the 
region covered by the wave becoming entangled with the outside region, see Fig. 1. This was 
dubbed an “entanglement tsunami” [3] with ve interpreted as the velocity of the tsunami 
wave front. It was further observed in [3] that for holographic models after local equilibration 
the normalized rate of growth 


^£(t) = 


1 dS^, 

'^eq-^S dt 


(1.5) 


^ See [4] for a proposed theory explaining this behavior. 
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appears to be bounded by (1.4) 


^s(^) < v_ 


,holo 
E 5 


t » £, 


eq • 


( 1 . 6 ) 


In the linear regime (1.1), 91s is constant given by ve, bnt in general it can have a compli¬ 
cated time dependence. Note that dSE./dt cannot be compared meaningfnlly across different 
systems or regions of different size as it generally scales with the geometric size of E and the 
number of degrees of freedom of a system. 91s was designed to provide an intrinsic measure 
for the rate of growth. With dimension of velocity, we expect that in relativistic systems 
91s should be constrained from cansality by some mnltiple of the speed of light. ^ 



FIG. 1. The growth in entanglement entropy can be visualized as occurring via an “entanglement 
tsunami” wave carrying entanglement inward from E. The region that has been covered by the 
wave (i.e. the orange region in the plot) is entangled with the region outside E, while the white 
region is not yet entangled. 

The simplicity and universality of the linear growth (1.1) begs for an underlying physical 
mechanism. In particular, it would be desirable to relate ve and 91s(t) to the speed of light, 
and to develop some intuition about the physical origin of the value of ve in (1.4). 

In this paper we first derive a formula which relates the tsunami velocity to the mntnal 
information of certain spacetime regions. The positivity of mntual information then leads to 
a proof that in relativistic theories ve is bounded by speed of light in all dimensions, i.e. 

ve<1, (1.7) 

^ Equation (1.6) implies < ^Seqds, which is reminiscent of the small incremental entangling theorem 
for spin systems [5] . See [6] for an attempt at a field theory approach to the problem. 
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although, as we will discuss in later sections, likely for d > 2 the inequality cannot be 
saturated. A different proof of (1.7) has been found by T. Hartman [7], which uses the 
monotonicity of relative entropy. 

We then consider various explicit models for entanglement propagation. Calabrese and 
Cardy [1] proposed a simple free particle streaming model to explain the linear behavior (1.1) 
in (1 + 1) dimensions. In this model, the injected energy density due to a global quench 
at t = 0 is assumed to create EPR pairs of entangled quasiparticles which subsequently 
propagate freely. At t = 0, entanglement correlations among quasiparticles are assumed to 
be local, ^ which eventually spread to large distances via free propagation of quasiparticles. 
In this model (1.3) comes from quasiparticles traveling at the speed of light. That (1.1) can 
be captured by such a simple model is remarkable and surprising. It appears to indicate that 
interactions do not play a role in the growth of entanglement, with the long-range entangle¬ 
ment of the final state solely coming from the spread of initial short-distance correlations. 
An indication that this success is likely an accident is that for more than one intervals, the 
model fails to reproduce the qualitative behavior of both holographic and CFT results [8-10]. 

The free streaming model can be generalized straightforwardly to a more general entan¬ 
glement pattern of quasiparticles and to higher dimensions. With no interactions, the wave 
function of the full system factorizes into those associated with each spatial point at t = 0, 
and the propagation of entanglement is determined by the entanglement measure /i[A] for 
a subsystem A in the Hilbert space of a single point at t = 0. In d = 2, due to the spe¬ 
cial kinematics of one spatial dimension, propagation of entanglement is independent of the 
choice of /u[A], and results from the EPR model of [1] are in fact general. This is no longer 
so in higher dimensions. Essentially all aspects of propagation depends on /r[A].^ 

A particularly interesting choice of ii[A\ is what we will refer to as a random pure state 
measure (RPS), for which the entanglement entropy for a subsystem A is proportional to 
its size, i.e. 

I^kps{A) = smin(I4, Vq), (1.8) 

where Va, Vq denote the volume for A and A (complement of A) respectively, and s is a 
constant (which depends on the specific system). The measure is motivated from the result 
^ At the onset of linear regime, the scale of entanglement correlations should be controlled by 1/T. It is a 

good approximation to treat them as strictly local when considering regions with 1/T. 

^ We study these aspects in detail for various choices of jii[A] in Appendix B. For example, while one again 

finds linear growth (1.1) at early times, the tsunami velocity ve now depends on /i[A]. 
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of [11] where it was found that the average entropy for a subsystem (with size smaller 
than half of the total system) in a random pure state is to a very good approximation 
given by its size. The rough intuition behind RPS is that for a sufficiently “equilibrated” 
system, all degrees of freedom are entangled with one another in an equal way, and thus the 
entanglement entropy of a subsystem is proportional to the number of degrees of freedom it 
contains. 

Using the strong subadditivity property one can show that the RPS measure in fact 
provides an upper bound for the propagation of entanglement among all free streaming 
models, leading to an upper bound for (thus an upper bound for ve) 

(1.9) 


where 


1 d = 2 


, ,free 


r(^) 



d = 3 


d = 4 
d = oo 


( 1 . 10 ) 


is calculated using the RPS measure. Note that is smaller than 1 for d > 2, because 
quasiparticles propagate in different directions. Comparing with (1.4) note that 


,free 


< V 


hole 


d> 3 


( 1 . 11 ) 


In other words, in higher dimensions, the spread of short-distance correlations limited by 
causality cannot account for the result (1.4) for strongly coupled systems. Thus interactions 
must play a role. The ratio decreases with d for d > 3, i.e. the higher the spacetime 

dimension, the more signiheant role of interaction is. In particular, as d ^ oo, —>■ 0, 

while ^ 1. 

This then motivates to introduce interactions in quasiparticle propagation. With inter¬ 
actions we then immediately face the problem of characterizing the quantum state of an 
interacting many-body system. Instead of confronting this very difficult problem directly, 
here we seek a qualitative understanding of how linear growth can arise in an interacting 
system and how interactions can enhance the spread of entanglement. For this purpose, we 
will consider the infinite scattering limit, i.e. we assume that scatterings among constituents 
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of a system are so efficient that the typical scattering time can be taken to be zero compared 
to the scales of entangling region size R and time t of interests. Holographic systems after 
a quench in the limit of equilibrium temperature T ^ oo can be considered as an example, 
as there the local scattering rate is controlled by T. In the infinite scattering limit, the 
evolution simplifies as interactions do not introduce additional scales into the problem. 

We propose a very simple model which applies the RPS measure (1.8) to certain spatial 
regions determined from the causal structure associated with the entangling surface. The 
use of RPS measure is natural in the infinite scattering limit as in this limit interactions are 
extremely efficient in redistributing entanglement. 

The model, to be referred to as the maximal RPS model, appears to capture gross features 
of entanglement spread, including the linear growth, of holographic systems in the T —>■ oo 
limit. In fact in d = 2, it does much better than expected.^ Applying it to a single interval 
in d = 2 we again find (1.3). For two intervals, the model precisely recovers holographic 
results. For three and four intervals, we find the model (and its slight variations) reproduces 
intricate entanglement patterns exhibited by holographic systems for a significant part of 
parameter space. For general d > 2, this model gives in all dimensions ve — 

We also show that the failure in reproducing holographic results for certain regions of 
parameter space for three and four intervals in d = 2, as well as = 1 for d > 2 can be 
attributed to the fact that our model can violate the strong subadditivity condition and/or 
does not take full account of causality constraints. 

Given the already remarkable success of the maximal RPS model, in an attempt to better 
understand how quantum information is organized in holographic systems, we construct an¬ 
other model for the evolution of entanglement entropy that is inspired by the tensor network 
description of state evolution after a quench. This interacting model satisfies many natural 
geometric criteria on the entropy function, including strong subadditivity, and reproduces 
the holographic results for multiple intervals in d = 2. However, it does not provide a 
stronger constraint than u® < 1 on the tsunami velocity for d > 2. It would be very inter¬ 
esting to find further criteria on the entropy function for the time evolution generated by a 
local Hamiltonian. 

The plan of the paper is as follows. In the next section we provide a proof that the 


® For multiple intervals, results from our model coincide with a phenomenological formula recently proposed 
in [9] to capture holographic results. 



tsunami velocity is bounded by the speed of light. In particular, we derive a formula which 
relates the tsunami velocity ve to the mutual information of certain spacetime regions. 
As a preparation for later discussions, in Sec. Ill we introduce the RPS measure for the 
entropy of subsystems and show that it is an upper bound for the entropy of systems of 
particles homogeneous in space. In Sec. IV we introduce free streaming models of entangle¬ 
ment propagation and discuss some explicit examples. A general upper bound on ballistic 
propagation of entanglement in these free streaming models is proven in Sec. V. The proof 
demonstrates that these models cannot account for the velocity of entanglement propagation 
in holographic theories. In Sec. VI we present an interacting model and apply it to various 
examples. It indicates that interactions can increase the tsunami velocity in higher dimen¬ 
sions and gives an upper bound for this velocity in relativistic theories. We also discuss 
shortcomings of the model. In Sec. VII we develop yet another description of entanglement 
evolution that is inspired by tensor network constructions. This model has the advantage 
of satisfying the strong subadditivity constraints and reproducing the right entanglement 
pattern for holographic theories in d = 2. In this context we also show the holographic result 
gives an absolute upper bound on the entropy after a quench for relativistic theories in d = 2. 
In Sec. VIII we further examine a relation discovered in Sec. II between tsunami velocity 
and mutual information in the context of various models discussed in earlier sections. Some 
more technical results and proofs are given in the appendices. In particular in Appendix B 
we work out many aspects of ballistic propagation of entanglement in general dimensions, 
including mutual information and finite volume effects, for different entanglement patterns 
of the initial state. 


II. PROOF OF A GENERAL UPPER BOUND ON THE TSUNAMI VELOCITY 

For d > 2, so far the linear growth (1.1) has only been established for holographic 
systems [2, 3]. In particular, both the linear growth (1.1) and tsunami velocity ve are 
independent of the shape of entangling surface S [3]. Assuming this behavior persists for 
general systems, here we prove that that the tsunami velocity is bounded by the speed 
of light. 

Since (1.1) is shape-independent, it is enough to consider S a straight hyperplane, i.e. the 
entangling region is a half space. Let us consider at time t a half-space W(t) whose boundary 
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is perpendicular to the a:i-direction and another half space hh(t + 5t) at time t + 5t. Since 
entropy in the quench evolution is translation invariant, we can take W{t+5t) infinitesimally 
displaced in the Xi coordinate such that the boundaries of these two regions W(t) and 
W{t + 5t) are connected by a strip of a null plane X (see Fig. 2). Consider the mutual 
information between W (t) and X 


I{W{t),X) = SiWit)) + S{X) - S{Wit) U X) 

= S{X) - + 6t)) - 5(fF(t))] 

= S{X) - VESeq^T^St , (2.1) 


where in the second equality we used S'(lF(t) U X) = S{W{t + 5t)) and in the third line the 
linear regime formula (1.1) for S'(VF(t)) and S{W{t + 5t)). 

Here we are interested in an excited state with vacuum entanglement subtracted, thus 
all quantities in (2.1) should be understood as so. 

Equation (2.1) then gives 


ve 


SiX)-I{W{t),X) 


( 2 . 2 ) 


and from non-negativity of the mutual information 

S{X) 

Ve < . • (2.3) 

It should be emphasized that here we are saying the mutual information with the vacuum 
part subtracted should also be non-negative. This can be justified since we are considering 
only contributions linear in Seq. We can take Seq large so the term proportional to Sgq in the 
mutual information has to be non-negative by itself.® 


W{t + 5t) 


W{t) 



FIG. 2. A half space at time t and another at time t+6t. The null region X connects the boundaries 
of these two regions. The horizontal direction is xi and vertical direction is time. Directions parallel 
to the boundary are not shown. 

^ Note that in the vacuum, there are divergences associated with the sharp corner between W{t) and X, 
which could ruin an inequality like (2.3). In the same spirit, peculiarities for the entropy of sharp null 
surfaces in interacting theories [12] do not apply here since in the present context all geometries can be 
thought of as having a minimum curvature radiuiOl/T. 
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FIG. 3. The small strip regions A, B, and C have all the same width. The null regions A' and B' 
are included in the causal development of regions A and B respectively. 

Now consider the three strip like regions A, B and C of figure 3. We take these regions 
to have small width 1/T <C fix = 25t <C t. Hence, from (1.2) the entropies for these regions 
are given by 

S{A) = S{B) = ^((7) = seq^sfix, (2.4) 

with A^ the area of the plane bounding the strips. From this volume law and the corre¬ 
sponding one for the strip AB, S{A[J B) = 2seq^sfix it follows that 

I{A,B) = S{A) + S{B)-S{AUB)^0. (2.5) 

The null regions A' and B' are included in the causal domain of dependence of A and B 
respectively. Then, because of the monotonicity of mutual information, the mutual infor¬ 
mation between the null strips A' and B' must also vanish 

0 < /(W, B') = SiA') + 5(5') - 5(H' U 5') < I{A, B) = 0. (2.6) 

Now given that 5(H' U 5') = S{C) and that by symmetry the entropy in the two null strips 
is the same, from (2.6) we find that 

St 

5(W) = 5(5')=SeqAsy = Seq^sfit ^ 5(X) = Seq^^fit, (2.7) 

i.e. the entropy of a small null strip is equal to the thermal entropy of its projection to a 
constant time slice. Plugging (2.7) into (2.3) we then get 


xe < 1, 


( 2 . 8 ) 


11 




which proves that for any relativistic system ve must be bounded by the speed of light. A 
different proof of (2.8) has been found by T. Hartman [7], which uses the monotonicity of 
relative entropy. 

Plugging the expression for S{X) into (2.2), we also find 




Equation (2.9) is an instructive formula which relates the deviation of a tsunami velocity 
from the speed of light to the entanglement between W{t) and the null surface X. In 
particular, it says for ve to be equal to 1, W(t) and X have to be unentangled. 


III. RANDOM PURE STATE MEASURE 

For the rest of the paper we investigate the propagation of entanglement in various free 
streaming and interacting models. Before doing that, we digress here to discuss the random 
pure state measure, which will play a crucial role in our later discussions. 

Consider a gas of particles on a spatial manifold M. in some pure state. We assume 
that the state is homogeneous, i.e. the particles are uniformly distributed on Xi. Now 
consider the entanglement entropy *S'(A) for a subregion A. We assume that S{A) satisfies 
the condition that as the size of A goes to zero 

limA(A) = sV:4, (3.1) 

A^O 

where s is a constant. We will comment on the motivations for this condition a bit later. 
We first show that given (3.1) the entanglement entropy *S'(A) is bounded by 

S{A) < )Urps(A) = smm{VA,VA), (3.2) 

where Va, Vq denote the volume for A and A (complement of A) respectively. We will refer 
to /rRps[A] as a random pure state measure. We apply the strong subadditivity condition [13] 
to regions A, B, C shown in Fig. 4 to get the inequality: 

A(A) + A(B U C) > A(A U C) + A(B). (3.3) 

Using that B, C are infinitesimal, we then have 

SiA U C) - SiA) < SiB U C) - SiB) = sVc , (3.4) 
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where we have used (3.1) on right hand side of the equality. This inequality holds for any 
region A and any other infinitesimal region C outside A. It implies that as we increase 
the size of a region the entropy cannot increase faster than the volume of the region times 
s7 Therefore, for any region A with Va <V/2 where V is the total volume we necessarily 
have S(A) < jURPs(A) = sVa, while for a region with Va > V/2 we have ^'(A) = S{A) < 
/^RPs(^) = /^RPs(^); proving our assertion (3.2). 



FIG. 4. Regions A,B,C on A4. 


Let us now elaborate on the condition (3.1). Suppose that as VT —?• 0, S'/! is instead given 
by Sa ~ with a 7 ^ 1. Then for Vg, Vc ~ e ^ 0, we find that 

SiA U C) - SiA) < SiB U C) - SiB) oc e“ . (3.5) 

For the case of super-volume behavior, a > 1, since e"/e —> 0 as e —>■ 0, we conclude that 
S'( 74 ) is identically zero for all A. Clearly, this is unphysical. 

For a < 1 , for small A, Sa can increase as Va which is faster than (3.1). This is 
indeed possible, but if this power law behavior continues to large sizes no volume law or 
equilibrium entropy after a quench can be achieved. In fact, SaIVa —t 0 for large volume 
if a < 1. Hence, our condition (3.1) simply means that for the quenched systems we are 

Similarly, from S{A) = S{A) it also follows that S'(^) — U C) < sVc, i.e. as we increase the size of a 
region the entropy cannot decrease faster than the volume times s. 
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interested in, which have large energy density with respect to subsystems sizes and times 
involved in the dynamics, we are assuming a local equilibrium has already taken place for 
small sizes and times of order 1/T. 

The condition (3.1) and the random pure state measure (3.2) are also reminiscent of [11] 
where it was found that the average entropy for a subsystem (smaller than half size of the 
total system) in a random pure state is to a very good approximation given by its size. 


IV. FREE STREAMING MODELS IN GENERAL DIMENSIONS 

A. Setup 

Here we consider a model of free propagation of entanglement after a quench at t = 0. 
We assume the system is homogeneous, isotropic and is in a pure state after the quench. 
We take the initial state at t = 0_ as unentangled, and require that in equilibrium (defined 
as t ^ oo) the entanglement entropy is proportional to the volume of a region. We make 
the following assumptions on the quench and the subsequent propagation of entanglement: 

1. The quench generates a large amount of short-distance entanglement correlations 
which subsequently propagate freely. We will assume that the quench time interval is 
negligible and the initial entanglement correlations can be taken to be local. In other 
words, at t = 0, each point acts as an independent source of entanglement correla¬ 
tions, which then spread freely, limited only by causality. At time t, the entanglement 
relations from T = 0 spread at most to the sphere \x\ = t. 

2. For simplicity, we will assume that the correlations are concentrated on the light 
cone. It should be straightforward to generalize the discussion to include entanglement 
correlations inside the light cone, which we expect not to change qualitatively the 
physical picture and the upper bound on the propagation derived below. 

3. There is no interactions/interference among light cones from different points. This 
implies that throughout the evolution different light cones can be associated with 
independent wave function factors inherited from those at the origin of each light cone 
at t = 0. Furthermore, entanglement correlations on each light cone do not evolve 
with time. More explicitly, consider a region A on a light cone from x = 0 with fixed 
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angular extension with respect to the origin. Denote the entanglement entropy of A 
with respect to the rest of the light cone as /i[A]. We assume that /i[A] is independent 
of time. 


Note that the model does not assume propagation of particles, only free propagation of 
entanglement correlations. The model is fully specified by the entanglement measure /i[A] 
on a light cone. When considering certain specific realizations of n[A\ it will often be useful 
to consider quasiparticles such as the EPR pair and GHZ block examples discussed below. 
We could also simply postulate a ii[A\ which satisfies all the properties of entanglement 
entropy as in the example of random pure state measure discussed below. 

Given that each light cone is independent, the time evolution of the entanglement entropy 
of the region A enclosed by a surface S can be obtained by summing over the entanglement 
entropy of the parts of the light cones intersecting this region, i.e. 

= J d‘^~^x iJ,[Lj^{x;t)], ( 4 . 1 ) 

where Z/s(t t) denotes the region(s) of the light cone with center x lying inside E at time t. 
fi[L] is zero if L is an empty set (i.e. no intersection). 

The entanglement measure ia[A] for a region ^4 on a light cone from x can be interpreted 
as the entanglement entropy for A by tracing out degrees of freedom outside A within the 
Hilbert space of x. It should satisfy all the properties of entanglement entropy for a pure 
state, including for example, 

^[A] = ^[A ], ( 4 . 2 ) 

where A denotes the complement of A on the light cone, and the strong subadditivity 
condition 

/i[A] + /i[ 5 ] > iJ,[A n 5 ] + iJ,[A U B] ( 4 . 3 ) 


for any region A and B. 

We also assume that measure fi[A] for small A (even with several disconnected pieces) 
is proportional to the normalized area of that region, which we already motivated 
around (3.5) in the last section. More precisely, for a region A included in a spherical 
cap of angular size A9 


lim 

A6»^>0 




= S, 




<^A 

^d-2 


(4.4) 
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Here s is a constant and ua denotes the volume of region H on a unit sphere and uJd -2 is 
the volume of a unit (d — 2)-sphere. In our current context, we will see in Sec. IV C 1 that 
imposing (4.4) is equivalent to the requirement of a final equilibrium state with a volume 
law entropy distribution. In fact, the final equilibrium entropy density is precisely given by 
the constant in (4.4), i.e. Seq = s. 

By definition, /i[H] and thus s has the dimension of an entropy density, i.e. 1/volume. 
Given that s is the only scale of the system, for a scalable surface E, on dimensional grounds 
we can write S's in a scaling form 


S^{t) = sR^-^f{tlR), (4.5) 

where i? is a characteristic size of S, or equivalently 

^<S's(t), (4.6) 

where AE is E rescaled by a factor A. This scaling relation is also satisfied by holographic 
systems in the large size and long time limit, as we will discuss more in Sec. VIB. Equa¬ 
tion (4.6) implies that for t small, when Sj: should be proportional to the area of E, Sj: must 
grow linear with t. At large times as t —?• oo, if the system has an equilibrium, i.e. f{t/R) 
has a well defined t ^ oo limit, then s must be proportional the equilibrium entropy density. 


B. Some examples 

It is instructive to look at some specific examples of /i[A]. 


1. Entanglement earned by EPR pairs 

One assumes that the quench creates a uniform density of independent EPR pairs of 
quasiparticles which are entangled within each pair and subsequently travel in opposite 
directions at the speed of light. The distribution of the directions of pairs is isotropic. This 
is a higher dimensional generalization of the model of [1]. Under these assumptions, at time 
t, a point X is entangled with another point y if and only if 

\x — y\ = 2t . (4.7) 
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If A consists of a region included in one half of the light cone it immediately follows that 


/^epr[^] — , (4.8) 

where s is a normalization constant and is the area of region A normalized by the area 
of the whole light cone (4.4). The normalization constant s can be written more explicitly 
in terms of the particle density n as 

s = niy , (4.9) 

where p the entanglement entropy within the pair from tracing out one of the particles. 
For general A, not necessarily included in a half light cone sphere the measure is more 
complicated and depends on relative locations of different parts of ^4. It can be written 
formally as 

= ^ f + n'), (4.10) 

l^d-2 JA JA 

where ft and n' denote unit vectors on a unit sphere with and Vt' the respective associated 
angular measure. We can give a more compact expression for this measure as follows. We 
define A' as the set of antipodal points in A, that is, A' is the set of unit vectors n such that 
—n G A. Then, as only those quasiparticles in A contribute to the entropy, whose pair (at 
the antipodal point) is in A, we have 

/^epr[^] = ■sCnnu' • (4.11) 

If A is included in one half of the light cone. An A' = A^ and we get back (4.8). 


2. 2m particles forming a GHZ block 

The EPR pair example has only bipartite entanglement among particles. We now consider 
an example with multipartite entanglement. One again assumes that the quench creates a 
uniform density of quasiparticles, but now particles from each point separate into uncorre¬ 
lated blocks each of which consists of 2m particles, with m an integer. Within each block, 
the 2m particles are entangled. To satisfy momentum conservation, for simplicity we will 
take that the 2m particles within a block come in pairs with back to back momenta. 

We consider the simplest entanglement relation that when tracing out any subset of 
particles within the block of 2m particles, one gets the same entanglement entropy p. This 
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( 412 ) 


is motivated by the so-called GHZ state for k qubits 

= ^ (| 0 )®‘ + | 1 )®‘) , 

which has this property and thus the name for this example.® 

Consider a region A which is included in half of the light cone, then we find that 

hGHzl^d; m] = — [1 — (1 — , (4-13) 

where s is again given by (4.9). To see (4.13), we note that : (i) ^ is the number of GHZ 
blocks at each point; (ii) from pairwise momentum conservation there can be at most m 
particles lying in A; (iii) as far as there are particles lying in A, the entanglement entropy 
for any particular configuration is always p; (iv) (1 — 2^^)"* is the probability that none of 
the 2m particles lies within A. 

For a general region, we have to calculate the probability that none or all particles are 
inside A. These probabilities are given by (1 — Cunu')"^ CunU' respectively. Hence, for 
any region A we have 

k'GHz[A; m] = — [1 - (1 - ^auA')”* - • (4-14) 

This formula gives back (4.13) for a region included in half of the fight cone, as 
and A n A' = 0 for this case. The case m = 1 reproduces the EPR result (4.11) since 
(Cuuu' — Cunu') = 2C^nA'- The formula (4.14) also obeys (4.2), which is most easily seen by 
rewriting it using 1 — ^^*4 A\J A' = Ar\ A' as 

hGHz[24;m] = — [l - ■ (4-15) 

3. Random pure state measure 

Another measure is the random pure state measure we discussed in Sec. HI, i.e. 

hRPs[-4] = smin(^^,^q) (4.16) 

for any region A and 5 is a normalization constant. While /iRps coincides with juepr when 
A consists of a region included in the half sphere, in general (4.11) is clearly different 

^ We note that because the tripartite information I 3 > 0 for GHZ states and holographic mutual information 
is monogamous [14], the GHZ pattern of entanglement cannot be realized using holographic geometry. 
Nevertheless, it is a simple example that holographic results can be compared to. 
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from (4.16). Note that (4.16) only depends on the area of a region on the unit sphere (or its 
complement) for any A. This is not so for both /liepr and jughz- As we showed in Sec. Ill, 
given (4.4), the RPS measure provides an upper bound for other measures. This will enable 
to us to establish an upper bound on 91s(t) and ve in Sec. IV. 

C. Evolution of the entanglement entropy 

We now use (4.1) to derive some general results for the evolution of ^'^(t). Here we will 
focus on universal features which do not depend on the specific form of /i[A]. In Appendix B 
we study many other aspects including mutual information and finite volume effects based 
on the specific examples of the last subsection. 

We will consider a compact surface S with a characteristic size R. In other words, R 
collectively denotes the curvature radii of S. 

1. Equilibrium value 



FIG. 5. Analysis of the late time behavior. The collection of points whose light cones 

intersect with E is given by a shell centered around E (green region). The approximate radius of 
the shell is t. For a fixed rj, the intersections of light cones with E from different values of y provide 
a foliation of region A enclosed by E. 
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The equilibrium value may be defined as 


= Sj:{t = +oo) . (4.17) 

As t —>■ oo, the size of light cones become much greater than that of S, i.e. t ^ R. Thus 
only small fractions of a light cone can be inside E. The collection of points whose light 
cones intersect with E will be denoted as As (t) and is given by a shell centered around the 
region, see Fig. 5. The integration over A/s(^) can be written schematically as 

j j dy, (4.18) 

where ff denotes directions tangent to shell, while y denotes the integration along the width 
of the shell as indicated in Fig. 5. Clearly, the precise shape of the shell will depend on the 
shape of E and for an irregular shell there may not be a preferred splitting in (4.18), but as 
we will see, such details are not important. 

Now fix an fj and consider the integral (4.1) over y. As we vary over the range of y, 
the corresponding provides a foliation of region A enclosed by E, see Fig. 5. In 

particular, for t —>■ oo, ^£( 77 , y] t) corresponds to a tiny part of the light cone from x = {ff, y) 
and we can approximate 

/«[As(??, 2 /;t)] = s^{LY:{ff,y;t )), (4.19) 

where we used (4.4) for infinitesimal regions. 

The normalized area for Lj:(ff,y;t) can be further written as 

= (4.20) 

where a{y) is area of Lj:(fi,y;t). We thus find the integral of (4.1) over y gives 

I dy fi[Lj:{ff,y;t)] = -/ dya{y) = , (4.21) 

where Vs = / dya(y) is the volume of the region A enclosed by E. Note that the above 
integral is independent of ff. Now integrating over ff, to leading order in large t approximation 
we simply get the area of the cross section of the shell, which is in turn given by the area 
of a sphere of radius t. Such a factor precisely cancels the denominator of (4.21) and we 
conclude that 

<S's(t —>• 00 ) = SeqVs (4.22) 
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with the equilibrium entropy density given by 


Seq = S . (4.23) 

For the RPS example s is simply a normalization constant and nothing more can be said. 
For the EPR and GHZ examples, the above equation can be further written in terms of the 
particle density n as 

Seq = nv , (4.24) 

which has a simple physical interpretation. Recall that for both examples v is the entangle¬ 
ment entropy for a single particle when tracing out the others. As t ^ oo, the entangled 
particles are separated by infinite distances and thus the entanglement entropy for any finite 
region is given by the particle number density times p. For a generic interacting system the 
equilibrium entropy density Seq is expected to coincide with the thermal entropy density, 
and (4.24) is also natural from that perspective. 

It is not difiicult to realize from this proof that a uniform volume law for infinitesimal 
regions of any shape on the sphere as in eq. (4.4) is also necessary to get the volume law at 
late times (4.22). 

2. Early linear growth 


Now let us consider the early growth, i.e. for t -C R. At such times, the radius of a 
light cone is much smaller than the curvature radius at any point of S, see Fig. 6. We can 
then locally approximate E as a straight hyperplane, with translational symmetries along 
directions tangent to E. The integrations in (4.1) can then be factorized into an integral 
along E, which simply gives a factor As (area of E), and the relative location y of centers of 
light cones with respect to E in the perpendicular direction. Then the early growth of the 
entropy will be determined completely by the measure p, applied to spherical caps. Let us 
introduce the notation 

hcap(^u) = h[A] (4.25) 

for the measure for spherical caps A as a function of their normalized area (note the area 
of a spherical cap determines it uniquely). Using this definition (4.1) can be written more 
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FIG. 6. Analysis of the early behavior. Left: For early times the curvature of S is much grater 
than the radius of the light cone, and S can be approximated by a hyperplane, shown on the right. 
Right: Variables used in (4.26) to calculate the early time evolution of entanglement entropy. 


explicitly as 

dyHct,Ai{ylt))=VESe^A^t, t<^R, (4.26) 

Jo 

where the factor of 2 comes from the domain of integration y G (—t, 0), 

2 

ve = — / dx ficapiCi^)), (4.27) 

Jo 

and ^{x) is the normalized area of a spherical cap for a unit sphere with angular spread 
defined by the perpendicular distance x = y/t from the center (see Fig. 6). In (4.26) we 
used the time independence of the entropy measure on the light-cone. In (4.27) we chose to 
normalize the quantity by the equilibrium entropy density Seq- 

Since (4.27) only involves a region smaller than half the sphere, Pepr and /Urps give the 
same result 


.EPR 


= V 


RPS 


= 2 


dx^(x) 


(^d-2 



dO sin'^-^ 6 = 


<^d-3 

~ dA2 

^d-2 — 



(4.28) 

For /tghz, a closed formula is not available, but one can easily obtain the numerical values 
for different values of m and spacetime dimensions, see Fig. 7. It is clear from the figure 
that for m > 1 and d> 2 


.GHZ 


< V 


EPR 


(4.29) 
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In Sec. IV we provide an upper bound for the speed ve for any entanglement measure. 



D. Quadratic growth before local thermalization 

In [3], it was found that for fast quenches in holographic theories, there is a period of 
quadratic growth before the linear growth, which sets in only after the local equilibration 
time teq- The local equilibration time is defined as the time scale by which local thermo¬ 
dynamics already applies with a thermal entropy density Sth, but long range correlations 
in which we are interested have not been established yet. For strongly interacting systems, 
various holographic studies [15, 16] indicate that the local equilibration ~ 1/T where T is 
the final equilibrium temperature, and by this time the entanglement correlations from the 
quench will have at most spread to a length scale 4q ~ l/T. For regions with size satisfying 
RT ;$> 1, we can treat tgq and 4q as zero, which is essentially what we have being doing so 
far. In other words, our above discussion should be interpreted as applying only after teq- 
More explicitly, in our setup we have assumed that at t = 0 entanglement measure fi[A] on 
light cones have already been fully established. We believe this assumption is reasonable 
only after local equilibrium has been established, i.e. after the quench has finished, it still 
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takes some time for a system to build up the local entanglement measure and that 

time scale can be interpreted to be 

We will now show that with some very simple assumptions, one can easily obtain quadratic 
growth of entanglement entropy with time. For definiteness of the discussion, we will use 
the EPR model as an example, although the discussion can be easily adapted for a generic 
measure n[A]. Consider equation (4.8), except that now the prefactor s is taken to be a 
function of time for t < tgq- We will take the simplest possibility: a linear function, i.e. 


s{t) 


, P t/ \ Lqq , 

f-eq ^ 

S t ^ tec[ • 


(4.30) 


For t > tgq the discussion is essentially the same as before, recovering the linear growth. For 
t < teq, the integral (4.26) becomes 

rt pt-to / 11 \ 7)EPR 

5's = 2Aj:— / dtoto / dy^ i -—— | = — Aj^t^ . (4.31) 

teq Jo Jo — to / ^teq 

In contrast, for holographic systems one finds [3] 


Ss = . (-i-32) 

where e is the energy density. This is not that different from (4.31) considering t^q ~ 1/T 
and that for a CFT e = s T. 


V. UPPER BOUND ON THE BALLISTIC PROPAGATION OF ENTANGLE¬ 
MENT 

In Sec. Ill we showed that the RPS measure (4.16) is an upper bound for all measures, 
an immediate consequence of which is an upper bound for the entropy of any region at 
any moment of time in models with ballistic propagation by an explicit geometric function. 
From (4.1) 

«Ss(t) < = Seq J d'^-^x min(CLs(«,t),Cz^) • (5-1) 

In this section we use this to prove an upper bound on 91s (^). 

First let us look at ve, which only involves spherical caps. Using 

/^cap(C(a^)) < /^Sp®(C(a:^)) = SeqC(a^), (5-2) 
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in (4.27) we find 


ve < Wg®® = 2 / dx ^(x) = 

Jo 


rf—' 

7h% 


(5.3) 


Note that the EPR measure (4.28) saturates the bound. We note that there is an alternative 
proof of 

hcap(^(^)) ^ •5eq^(^) ) (^-4) 


without using (3.2). In fact one can use the strong subadditivity condition to prove a 
stronger inequality 

f‘Lp(O<0, (5,5) 


i.e. /Ucap(0 is a concave function. We give the proof of (5.5) in Appendix A. Because a 
concave function always lies below any of its tangents, and j^cap(O) = ^'eq, equation (5.4) 
then follows. 

We now show that the normalized rate of growth of entropy 93s (^) is bounded by at 
all times. From (5.1) we get 


93s (t) 


1 dAs(t) ^ f ^d-l^ 

SeqAs dt As J dt 


(5.6) 


Therefore a bound on normalized grow rate 


93s (t) < 


(5.7) 


would follow from a purely geometric inequality 


d 


d-i 


X 


d^L^{x,t) 


dt 


< 


( 6 . 8 ) 


We now prove this last inequality (5.8). Let us write an integral representation for the 
normalized area: 

[ d'^~^y^i\y\ -l)Oi:ity + x), (5.9) 

^d-2 J 

where 0s (T) is the characteristic function of S, which takes the value 1 for z inside E and 
0 for T outside of it. The time derivative is given by 

d^L^ix,t) ^ f f daX^{y-n{a))5{ty + x-x{a)), (5.10) 

dt UJd-2 J Js 
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where a denote collectively a set of coordinates on the entangling surface S, n{a) is the unit 
vector normal to S, and x{a) the position vector on the surface. Then we have 


/ 







dt 



1 

[ d^-^y5{\y\-l) 

/ daX^ {y-n{a))S{ty + x 

^d-2 . 

1 

ds 

— 1 

‘^d-2 J 


1 daXm \y-n{a)\ 5{ty + x 


x{a)) 
x{a )), 


(5.11) 


where we took the absolute value inside the integrals. By performing the integral over x, we 
can get rid of one of the delta functions, and we obtain 


d 


d-l 


X 


{x,t) 


dt 


< 


—f X ^yS{\y\-l) f daX^\y-n{a) 
^d-2 J Je 


(5.12) 


By using the rotational symmetry of the integral over y in (5.12), we can make the replace¬ 
ment n{a) n, with a fixed (and arbitrary) unit vector n, thus the integral over a can be 
evaluated to give 


I 


d'^-^x 


(x,t) 

dt 


< f—*— [ X ^y5{\y\ - 1) \y-n 

\^d-2 J 


(5.13) 


where in the last line we recognized (4.28). This completes the proof. 

Note that the bound for the growth rate can only be saturated if both inequalities (5.6) 
and (5.8) are saturated. The inequality (5.6) is only saturated for all light cones, if: (i) 
the measure is equivalent to the RPS measure; (ii) all light cones whose intersection area 
with the region A bounded by E is less than half the sphere volume are increasing their 
intersection with time; (iii) all light cones with intersection greater than half the sphere 
volume are decreasing their intersection with time (the latter two conditions are satisfied by 
any shape at t = 0). Even for the RPS measure this is not the case for long enough times 
(except for a planar entangling surface). We also note that 9d£(t) can be negative in some 
circumstances, as shown in Appendix Bid. 
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VI. INTERACTING MODELS 


In this section we present an interacting model.® Although the model of [1] captures 
the time evolution of entanglement entropy for a single interval precisely in two dimensions, 
we saw that in higher dimensions, and qualitative differences also arise for 

entanglement entropy for multiple intervals in two dimensions [8, 9]. Furthermore, in free 
propagation models the spread of entanglement depends sensitively on the entanglement 
pattern of the initial state as we have studied in detail in Sec. IV and Appendix B. In 
contrast, in holographic systems, as emphasized in [3], the linear growth emerges after a 
system has locally equilibrated with thermodynamical concepts such as temperature and 
entropy density already applying at scales of order 1/T. This implies that by the time the 
linear growth emerges, the details of the initial state should have already been largely erased 
by interactions, and the linear growth must be a consequence of interactions. 

An immediate consequence of interactions is that the quantum state of the system can 
no longer be described as a tensor product of those resulting from each point at t = 0. As 
a result, our fundamental equation (4.1) breaks down, and one has to face the problem of 
characterizing the quantum state of an interacting many-body system. Instead of confronting 
this very difficult problem directly, here we seek a qualitative understanding of how linear 
growth can arise in an interacting system and how interactions can enhance the spread of 
entanglement. We will present a very simple model, which we believe gives an upper bound 
on the spread of entanglement in interacting theories. In particular, it gives u® = 1 in any 
dimensions saturating the bound (2.8) in Sec. II. 

A. The effect of a scattering event 

To develop some intuition on the possible consequences of interactions on the propagation 
of entanglement, let us first consider the effects of a single scattering in the EPR model of 
quasiparticle propagation in (1 + 1) dimensions. Again we consider a single interval A of 
length 2R, with various different situations of scattering depicted in Fig. 8. Diagrams (a), 
(b), (c) can happen for t < f, (a), (b), (c), (d) for R > t > f, and (a), (b), (d), (e) for 
^ Interactions between quasiparticles and their effect on the tsunami velocity was recently investigated in 
(1 + 1) dimensions by Cardy in [17]. The interactions he considered break the scale invariance of the 
theory, whereas our setup preserves it. 
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t > R. We will approximate each particle as a qubit and treat the scattering as a unitary 
transformation. In other words, after scattering, particles 2 and 3 maintain their original 
directions but the resulting state is related to the product state before the scattering by a 
unitary transformation. Note that an arbitrary unitary can encode not just forward, but 
backscattering as well. The latter is implemented by a unitary that swaps the wave functions 
of 2 and 3. 

More explicitly, with the notation 

|1) = |00)23, |2) = |01)23, |3) = |10)23, |4) = |11)23 , (6.1) 

the scattering of 2 and 3 can be described by 

K) ^ 10 = U\i) = U^,\j), i,j = 1,2,3,4, (6.2) 

with Uij a unitary matrix. Without loss of generality the state for 1,2,3,4 before the 
scattering can be taken to be 

~ 1—rn—|2 (I®®)i 2 + cr|ll)i2) ® (100)34 + tt|ll)34) ( 6 - 3 ) 

1 + lap 

with a parametrizing the entanglement among the pair. The entanglement entropy obtained 
from tracing out one of the particles in a single pair is 

u = logp„ - (1 - Pec) log(l - Pec), Pa = , 10 • (6.4) 

1 + |a;| 

The state after the scattering is 

I?/;/) = (1 0 [/0 1)10) • (6.5) 

For different situations depicted in Fig. 8, we are interested in different reduced density 
matrices: 

(a) In this case the relevant quantity is S' 13 , i.e. the entropy for the reduced density matrix 
Pi 3 from tracing over particles 2 and 4. We will denote the entropies corresponding 
to (6.3) and (6.5) as and respectively; Si^ = 2z/. Clearly for |q;| = 1, i.e. when 
an EPR pair is maximally entangled, Os is always smaller than the maximal possible 
value S'jg = 2 log 2 for any U ^ 1. Physically, some of the initial entanglement between 
(1, 2) and (3,4) is now shared between 1 and 3, and 2 and 4. In the other extreme, for 
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FIG. 8. The effect of one scattering event on the entanglement of region A. The scattering is 
represented as the red dot, where a unitary operator U acts on the Hilbert space of particles 
(2,3). The labeling of the indistinguishable particles (2,3) after the scattering is arbitrary, but it 
is convenient to choose that the particles maintain their original directions. For a single interval 
A we show all the nontrivial scenarios, and we analyze the consequences of each configuration in 
the text. 
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\a\ <C 1 (or equivalently |q;| ^ 1 ), i.e. when the original pair is only slightly entangled, 
the scattering between 2 and 3 should generate new entanglement between the two 
particles and thus enhance the entanglement of (1,3) with the rest of the particles, 
i.e. we expect for generic U. For a in between, some U could enhance 

the entropy, and some could reduce it. The explicit expression for S'lg is given in 
Appendix C. 

(b) The relevant quantity is Si. Since U only acts on the complement of 1 , clearly Si is 
not modified by the scattering. 

(c) The relevant quantity is Si = S '123 and as in (b) the value is not modified. 

(d) The relevant quantity S '23 is again not modified, as U acts on the subspace of 2 and 3. 

(e) The relevant quantity is S^. When a pair is maximally entangled (|a;| = 1), S^ is not 

modified for any [/, as the entanglement between 3 and others is simply redistributed. 
For other values of a, the situation is a bit similar to that of S '13 discussed above in (a). 
For a 1 or tt ;$> 1, the entanglement will be enhanced for generic U. For generic 

values of a in between, then depending on U, the entropy can be either enhanced or 
reduced. The explicit expression for S^ is given in Appendix C. 

One should be able to use the above analysis to construct a dilute gas model to obtain 
quantitative results, as there are only a small number of scatterings. We will leave this for 
the future. Here we consider a model for the infinite scattering limit. 


B. The infinite scattering limit 

When including multiple scatterings, one needs to consider states associated with more 
and more particles and clearly we lose control very quickly. Nevertheless based on intuitions 
obtained from single-scattering results, we will see that one could still draw some qualitative 
conclusions about the general interacting case. 

On dimensional grounds we can write the evolution of entanglement entropy for a region 
A enclosed by S in a form 


Sji{t) = s,^R<^-^f{t/R,^/R,TR), 


( 6 . 6 ) 
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where T is the equilibrium temperature, Seq is the equilibrium entropy density, and ^ col¬ 
lectively denote other length scales in the system. In strongly interacting systems such as 
holographic theories, the local scattering rate is controlled by 1/T. For weakly interacting 
systems, the mean free path is typically controlled by other scale(s) ^ ^ 1/T. To simplify 
the analysis, we will work in the regime 

7-0, TR^oo, = finite, (6-7) 

and assume that the scaling function / (4.6) has a fini te limit in this case: 

= SeqR'^~'^fit/R) . ( 6 . 8 ) 

We refer to this regime (6.7) as the infinite scattering limit. For example, holographic 
systems after local equilibration are governed by it. Note that (6.8) is of the same scaling 
form as (4.5) valid in the free propagation case, even though the underlying physics is very 
different. In free case (4.5) is a consequence of no scattering and s is determined by the 
initial state, while in (6.8) is essentially a consequence of infinite number of scatterings (as t 
is infinite compared with any scattering time) with Seq determined by dynamics. As already 
mentioned below (4.6), the scaling form (6.8) implies that if there is a regime that Ss is 
proportional to the area of E, then the time dependence must be linear. 

We again assume that the quench generates a finite density of identical particles, which 
then subsequently propagate at the speed of light isotropically, and we allow an arbitrary 
number of scatterings. We will assume that on average scattering events are isotropic and 
homogeneous in space, implying that both incoming and outgoing particles are uniformly 
distributed in all directions. As in the (1 -|- l)-dimensional example of last section, scattering 
events will be treated as unitary transformations on all particles (which are assumed to be 
identical) that are at the same point at a given time. The labeling of particles after a 
scattering event is again arbitrary. Given the isotropy we can choose the labeling so that 
particles will not change direetions after scatterings. This means that we can trace the whole 
spacetime trajectory of a particle from its origin even in the presence of interactions. The 
ability of tracing a particle trajectory will play an important role in our discussion below. 

At any given time, to calculate the reduced density matrix for A (and thus the entan¬ 
glement entropy), we need to simultaneously trace over all particles lying outside A, rather 
than restricting to a subspace like in the non-interacting or single-scattering case. Now note 
the following: 


31 


1. The situation in (b) of last subsection can be immediately generalized to conclude that 

scatterings that happened in the past domain dependence of A (complement of A), 
i.e. in are not relevant, as they amount to unitary transformations in that 

do not change Pa-^^ In Fig. 9, where we depict the time evolution for one interval in 
(1 + 1) dimensions, these are regions shaded in red. 

2. Similarly as in (c) and (d), scatterings among particles which fall inside region A are 
also not relevant, as such scatterings act by unitary transformation on ?{a and thus 
will not change pa- In other words, scatterings in region T>_{A), shaded green in Fig. 9 
can be neglected. Thus particles that spend their whole life in the green region do 
not give rise to any entanglement. Note that for t > R, the green region no longer 
intersects with the t = 0 spatial manifold. 

3. Situations like (a) and (e) corresponds to scattering between particles one of which 
falls into A and the other falls into A. We will refer to these as effective scatterings. 
As in the single scattering case effective scatterings do affect entanglement. 

The above discussion shows that we only need to consider particles originated from the 
region 

Rf{t) =M- (V-(A) DM)- (V-(A) n M) , (6.9) 

where Ai denotes the full spatial manifold ai t — 0.^^ Furthermore, only those scatterings 
of these particles that take place in the white regions in Fig. 9 are relevant. This implies 
that disconnected regions of A/'(t) can be treated independently of one another. 

In the regime (6.7) all the particles will have scattered essentially an infinite number of 
times. In such a situation we expect that any memory of the initial state will be forgotten, 
and we can simply assign a geometric measure for the entanglement entropy. Since now the 
relevant Hilbert space is that for all the particles in A/", we will simply postulate a random 
pure state measure for the entanglement, i.e. 

SA{t) = I'eq [Afi{t )], Na [^(t)]) , (6.10) 

i 

NA[Afi{t)]= f nA{x,t), [^^(t)] = [ nA{x,t), (6.11) 

J Mi {t) J Mi {t) 

Recall that the past domain dependence V_ {A) of A is defined as the spacetime region where all future- 
extended causal curves pass through A. 

Note that the region A/e(^) we used in the free streaming model is in general a subset of A/’(t), see Sec. VIE 


for further discussion. 
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FIG. 9. Entanglement entropy of an interval Ain (1 + 1) dimensions in the interacting model. The 
top figure shows a time before saturation, while the bottom figure applies for times after saturation. 
The dashed purple lines show particle trajectories, and their intersections are scattering events. 
The green region is V-{A) and the red region is V-{A). As explained above scatterings that take 
place in these colored regions do not change the entanglement entropy. The point P at distance x 
from V-{A) and the left moving particle emanating from it will play a role in the discussions in 
Sec. VIE. 

where z/eq may be interpreted as average entropy per particle, and Na [Afi{t)] is the number 
of particles originated in Afi{t) that fall into A at time t. The number of such particles is 
given by the integral of the density of particles nA{xA) that originated from x and fall into 
A at T Taking the smaller value of the number of particles falling into A and A has the same 
rationale as the earlier postulate of random pure state, and ensures Sa = Sa- The sum i is 
over the disconnected components of A/", i.e. J\f = UiAfi. As mentioned earlier, disconnected 
components of J\fi should be treated independently and thus summed separately. 

At small t, the number of Afi is always the same as the number of disconnected boundaries 
of A. As time evolves, different Afi can join each other. Eventually all Affs will merge into 
a single connected region Af{t). At sufficiently late times, for A in an infinite space, it will 
always be the case that Na [Af{t)] < [Af{t)] as A is infinite. From homogeneity we can 
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then conclude that 


Sa = Seq^A) ^eq = VeaJ^i t Sufficiently large, (6.12) 

where n is the number density. 

We move on to analyze some examples. 


C. (1 + l)-dimensional examples 


Let us first consider various examples in (1 + 1) dimensions. In this case, the time 
evolution of entanglement entropy that we obtain (6.10) coincides with a phenomenological 
formula recently proposed to describe holographic results in [9]. In the limit R,t ^ 1/T 
the holographic result for N intervals can be obtained from a very simple minimization 
procedure: take all possible pairings of the left interval endpoints ■ ■ ■ ,^n} with the 

right interval endpoints {ri,..., rjv} and connect them by an extremal surface in the bulk. 
For the area of one extremal surface we get the one interval result 


ft {t < R) , 

5'interval(L-R) = 2Seq \ (6.13) 

fi? {t>R). 


We have to add up these contributions and minimize over the pairings: 


Snit) 


min 

a 


N 




interval 



(6.14) 


\ / _ 

where cr is a permutation. Note that for different times a different permutation may realize 
the minimum, which corresponds to the change of dominance of the extremal surfaces in 
holography. 


1. One interval 

Consider a single interval of length 2R in (1 + 1) dimensions, see Fig. 9. For t < R, we 
have two disconnected A/) each of width 2t, and for any point ua = Ua = ^, where n is the 
total particle density. We then find that 

Ti 

*^A(f') ^eq X 2 X — X 2t 25eqt, S^q = t R . (6.15) 
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For t > R, there is only one connected Rf, and 


Na [A/'(t)| = 

/ Tl 

/ nA{x, t) = 2 X — X 2R = 2nR , 

( 6 . 16 ) 


JAfit) 2 


Na ^(t)] = 

/ ttq(a:, t) = 2nR + 2{t — R)n = 2nt , 

( 6 . 17 ) 


JM(t) 



thus we conclude that 

SAit) = 2Seqi? = t> R. (6.18) 

We then find that 

ve = 1 . (6.19) 

The time evolution of entanglement is in agreement with (6.13). 

2. Multiple intervals 

Let us first explain how the calculation goes in an example of two intervals. For definite¬ 
ness, consider intervals of lengths 2Ri, 2 R 2 separated by a distance L with L < 2Ri < 2 R 2 , 
see Fig. 10. As in Fig. 9, green light cones indicate T>_{A) and red ones At a given 

time, indicated by the horizontal blue line we decompose Af{t) into disconnected pieces. At 
the time indicated in the plot we have two disconnected pieces. We then count the number 
of particles that end up in A and those that end up in A for each A/), and take the minimum 
of the two numbers. 

From Fig. 10 and elementary geometric considerations one readily sees that the number 
of particles from each A/) that end up in A and A is: 

Na [A/i] = n vol (A/) n A ), [A/)] = n vol (A/) fi A) , (6.20) 

and (6.10) can be rewritten in a simpler form: 

SA{t) = Seq min [vol n A ), vol fl A)] . (6.21) 

i 

Equation (6.21) was proposed recently in [9] to capture the time evolution of entanglement 
entropy in holographic systems, motivated from the picture of entanglement tsunami [3]. In 
this interpretation, A/'i(t)’s are given by the regions covered by the tsunamis,and one again 

Entanglement tsunamis originate from the boundaries of entangled regions and propagate in both direc¬ 
tions. 
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FIG. 10. Top: Explanation of how to calculate the entanglement entropy in the interacting model 
for two intervals of length 2Ri = 2 and 2i?2 = 3 separated by a distance L — 1. Light rays 
starting at the entangling surface partition J\fi into multiple pieces characterized by where the left 
and right movers end up. In the figure we label each such piece by two letters, with the left letter 
standing for left movers, the right one for right movers. For example AA means left movers end 
up in A while right movers in A. Bottom: Time evolution of entanglement entropy and mutual 
information for the same two intervals. The time slice considered in the top row is drawn by a 
solid blue line. The results agree with what one gets from the holographic result (6.14). (In the 
holographic calculations of e.g. [18] one sees the smoothed out versions of these plots, as they do 
not take the i?, L ^ l/T limit.) 

applies RPS to each Afi which amounts to taking the smaller volume between those of A 
and A regions within an Ni. As time evolves, different Afi regions join when their respective 
tsunamis meet. 

One can readily check that (6.10) (and equivalently (6.21)) reproduces holographic results 
for two intervals (6.14) for all parameters, and L. 

For three and four intervals, we find the model still reproduces precisely the rather in- 
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tricate holographic results for a significant part of the parameter space. With a random 
sampling of the parameter space with a large number of examples, for 3 intervals about 
21% of the examples and for 4 intervals 36% deviate from the holographic results.Even 
when (6.10) (and (6.21)) deviates from the holographic results, the overall trend is still quite 
similar, but it can give unphysical answers. We give an example in Fig. 11, where the en¬ 
tanglement entropy develops a discontinuous jump at some time. The mathematical reason 
for the jump is as follows. Let us denote the time of the jump as tc- At tc — e (with e ^ 0) 
we have two disjoint regions A/i and for A/j, [A/j] is smaller than A% [A/j], while for 
A/2, Na [A/2] is smaller, thus 

S{tc -e) = u^^iNA [AAi] + Na m) • (6.22) 

At tc + e, A/i and A/2 join into a single region and we should count the particles of Na and 
Na for A/i U A/2, i.e. 

A(% + e) = iy,^mm{NA[Afi] + iVq[A4], /Vn[AAi] + Na^) • (6.23) 

There is a discontinuity from (6.22) no matter which is chosen in (6.23). This phenomenon 
is likely generic as one increases the number of intervals: whenever two regions where A and 
A dominate respectively join together, there will be a discontinuous jump. In Fig. 11 we 
also plot the corresponding holographic result which is continuous and is smaller than (6.21) 
for a period after 4- We will further discuss the physical origin of the jump in Secs. VIE 
and VIF. 

Given its simplicity, it is quite remarkable that the maximal RPS model manages to 
reproduce intricate holographic results for multiple intervals for a significant part of the 
parameter space. As already alluded to in the Introduction and at the beginning of this 
section, the ballistic picture of the spread of entanglement discussed in Sec. IV-V of this 
paper cannot account for the holographic results for more than one interval in (1 -|- 1)- 
dimensions. This was analyzed in detail in [8, 9], and we do not repeat this comparison 
between the results of the ballistic and scattering pictures. Recently, a CFT analysis showed 
that rational CFTs behave according to the ballistic model, while non-rational CFTs are 

We sampled the parameter space by fixing the leftmost and rightmost boundary points and by throwing 
the other boundary points between them randomly using the uniform distribution. Our sample size was 
500. 
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FIG. 11. Top: Setup with three equal length intervals. We chose R = 1/2, Li = 1/20, L 2 = 3/2. 
We drew the past domain of dependence for our setup, with times before and after the jump 
marked by solid blue lines. Bottom: Time evolution of entanglement entropy for our setup with 
a jump between the blue solid lines. On the right hand side we illustrate how holographic CFTs 
behave. We drew the contribution from all possible locally extremal surfaces (i.e. the contribution 
from every possible permutation in (6.14)) by dotted green lines and the smallest one among them 
by a solid green line. Note that for different times surfaces connecting different endpoints of the 
intervals dominate. The result from (6.21) is again drawn by orange. 

expected to interpolate between the free streaming and the holographic behavior, which is 
reproduced by (6.21) [10].^^ 

Finally, note that there is a simple “phenomenological fix” to the discontinuity problem, 
as follows. When two RPS regions with different dominance join, say Afi with A and A /2 
with A as in the example of Fig. 11, immediately after the two regions start overlapping 
we still keep A// and A /2 as independent, i.e. the contribution from them is still given by 
5'rps(A/i) + aSrps(A/2) rather than the bigger value aSrps(A/i +A/ 2 ). As time evolves we merge 
them into a single RPS region when aSrps(A/i + A/ 2 ) = 5'rps(A/i) + aSrps(A/ 2 ). If before these 
The holographic behavior is universal in large central charge CFTs with a sparse low-lying operator 
spectrum. 
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two regions merge, other RPS regions start overlapping with either of them, we follow the 
same rule in deciding whether they should merge with other regions. Sampling over the 
parameter space, we found that in the improved model the failure rate in reproducing the 
holographic answer (6.14) is only 2% for three intervals and 8% for four intervals. 


D. Higher dimensions 

Let us first consider early times t R, for which we can approximate the boundary as 
a straight-line, then A/'(t) has is a strip of width 2t with E lying in the middle as indicated 
in Fig. 12. For this geometry the total number of particles from J\f falling into A or A are 
the same, so we can choose either of them. It then immediately follows from (6.10) that 

ST.{t) = X 2A^t = SeqA^t ve ^ I, (6.24) 

where in the above equation 2AY,t is the volume of Afit) and np^q is divided by two, as 
exactly half of all particles from Afit) will fall into region A because of isotropy. Thus, in 
this model the tsunami velocity ug is precisely given by the speed of light in all dimensions! 

It is interesting to contrast this computation with the earlier free propagation calculation 
of ve in Sec. IV C. We take the EPR (or equivalently the RPS) measure for the free streaming 
model. In both the free and interacting models we associate a measure to points on the t = 0 
time slice. Consider the contributions from points P and Q in Fig. 12. In the earlier free 
propagation calculation, for light cones originating from P or Q we took the area of the 
smaller spherical cap of the intersection of the light cone with S, while for (6.10), we simply 
use the spherical caps inside A, which for P is the bigger spherical cap. This thus leads 
to an enhancement in ve- Physically, in free particle models entangled particles originating 
from the same point and ending up in A do not contribute to the entanglement, while in 
the interacting model the initial entanglement pattern is forgotten and all particles ending 
up in A contribute. 

Now consider a spherical region of radius R. Region Af{t) is now the white annulus region 
indicated in Fig. 13. This annulus is cut into two parts by S, with the inner part having 
smaller volume. Then applying (6.21) gives that the entanglement entropy is equal to the 
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FIG. 12. Illustration of the computation of entanglement entropy for early times. We suppress the 
time direction and look at the configuration from “above”. Region A is to the left of the entangling 
surface S. The green region is D_(A), the red region is and the white region in between 

them is Points P and Q are on the t = 0 time slice and the purple light cones show where 

the particles that started out in P and Q end up at time t. 


volume of the inner white annulus in in Fig. 13: 


S{t) = Seqi?"-' 




VbcI-I 


('^ < 1 ), 
("^ > 1 ) > 


(6.25) 


where we have introduced r = t/R and denoted the volume of the unit ball by V^d-i = 

By taking the early time behavior we see that ve = I- This result can be contrasted with 
the ballistic propagation calculation for spheres presented in Appendix Bib. 


E. A family of RPS models 

Recall that in the free streaming model the wave function \'ijj) factorizes into those of each 
spatial point at t = 0, i.e. 

I'lp) = OxiV'x) (6.26) 

A nice check of (6.21) is that (6.10) gives the same result. 
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FIG. 13. A disk in d = 3 before saturation. The green region is P_(A), the red region is D_(A), and 
the white region in between them is Af{t). According to (6.10) by summing up the contributions of 
light cones like the one drawn with purple, we obtain the time evolution of entanglement entropy. 
Of course, it is a lot easier to use the simplification provided by (6.21). 

and the upper bound on the entanglement propagation is achieved when using RPS for 
each \^^). When including interactions clearly (6.26) does not apply. Nevertheless due to 
constraints from causality, within a finite interval t not all degrees of freedom can interact 
with one another. The basic idea behind (6.10) is that at time t, the full wave function can 
be factorized based on the casual structure of A, i.e. 

1^) = (®#M(t))) ® (• • •) > (6-27) 

where • • • denotes the factor of the wave function which is irrelevant for the entanglement 
of A. We then apply RPS to each 

Instead of (6.27) one can in principle consider a finer partition of Af{t) than connectivity, 

1^) = ® (•••) > (6-28) 
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where 


UaMa{t) = j^{t) = , (6.29) 

and for any a, there exists an i such that A4a{t) C We can obtain a general class of 

RPS models of entanglement propagation by applying RPS to each Aia, i-e. 

= z/eq min {Na [Wfa(t)], Na [Wfa(^)]) • (6.30) 

a 

The free streaming RPS model (6.26) is a special case of (6.28) with Aia given by a 
point.In other words, the free streaming RPS model is the hnest division of Af{t), it is 
the minimal RPS model. In contrast, the model (6.10) is the coarsest division, and thus we 
will refer to it as the maximal RPS model. 

Now using the same argument as demonstrating the discontinuity in (6.22)-(6.23), we 
can show that the entropy always decreases with a finer partition of Af{t). To see 

this, let us consider dividing some block Ai within a partition into Aii U Ai 2 = Ai. We 
then have 


Nm (A) = Nm, (A) + Nm, (A) , Nm (A) = Nm, (A) + Nm, (A) (6.31) 

and thus 

mm{NA4{A), NAi(A)) = min(A(A4i(^) A Nj^^iA), NAi^(A) + Nj^^iA)) 

> mm{NMAA), Nm,(A)) + mm{NM2{A), Nm2{A)) . (6.32) 

The fact that a subdivision decreases the entropy implies that with the class of all RPS 
models, the free streaming and the maximal RPS model provide respectively the lower and 
upper bounds, 

Sf%t) < . (6.33) 

Furthermore from the discussion of Sec. Ill, the RPS measure provides an upper bound for 
the entanglement entropy among all possible entanglement measures with a given partition. 
We thus conclude that the maximal RPS model (6.10) should provide an upper bound on 
entanglement propagation for all relativistic interacting systems after a global quench. This 
conjecture is consistent with our empirical observation that when the maximal RPS result 

Note that in free streaming model only the points in A/s(t) S Af{t) contribute, whose light cone can 
intersect with A. 
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deviates from the holographic one (6.14), it always lies above it.^^ That for d > 2 the 
maximal RPS gives ng = 1 is also consistent with the result of Sec. II. 

F. Further discussion of the maximal RPS model 

We expect that the RPS measure applies when all degrees of freedom in the relevant region 
(i.e. within each disconnected A/)) have fully “equilibrated,” i.e. have interacted sufficiently 
with one another. Otherwise it provides an overestimate. Consider at some time tc, there 
are two regions A/i and A /2 joining into a single connected region. (6.10) then dictates that 
at tc + e (with e —>■ 0) we must apply the RPS measure to the whole A/i U A^- But physically 
in going from tc — e to tc + e, there is just not enough time for this “total equilibration” 
to happen. When A/i and A /2 are dominated by A and A respectively before joining, such 
“lack of equilibration” will lead to a discontinuity, as in the example of Fig. 11. When A/ 1,2 
are both dominated by A (or A), or for one of them A and A give equal contributions, there 
will not be a discontinuity. Nevertheless, one may have expected that even in such cases 
the “lack of equilibration” may also lead to deviations from holographic results during the 
subsequent evolution after the joining. It is then rather curious that we do not observe such 
deviations at least in our sampling of the parameter space. This indicates that holographic 
systems “equilibrate” remarkably efficiently. 

The discontinuity in the example of Fig. 11 also means that the strong subadditivity 
(SSA) condition is violated, as SSA implies that the time evolution should be continuous as 
can be seen from the results of Sec. II. The violation of SSA can also seen from the behavior 
of the entanglement entropy at a fixed time as follows. Consider at time t = tc — e another 
three-interval region A which differs from A only by having a slightly smaller L 2 = L 2 — S 
with h ^ 0, e — 7 - 0, h > e. Then at time tc — e for A, the corresponding A/j and A /2 have 
already joined, hence at tc — e the entanglement entropies for A and A differ by a finite 
amount despite the fact that the two regions only differ infinitesimally, which violates SSA 
as expressed in (3.4). Thus, if a model can be constructed which guarantees SSA, then such 
discontinuities cannot arise. 

The issue of “lack of equilibration” for A/'(t) becomes more significant in higher dimen- 
In Sec. VIIC we prove that in (1 -|- 1) dimensions holographic theories give the fastest possible entan¬ 
glement spread. This then implies that whenever the maximal RPS result deviates from holography, it 
overestimates the entropy. 
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sions. Recall the calculation of ve illustrated in Fig. 12, for which N{t) has only one 
connected component. Clearly degrees of freedom far separated in directions parallel to 
the boundary of the region cannot be in direct causal contact. Phrased slightly differently, 
the maximal RPS model fails to take into account of causal constraints along longitudinal 
directions. Note that in this particular setup SSA does not appear to be violated. Thus 
longitudinal causality constraints (which arise only in d > 2) should be considered as an 
independent requirement from SSA. 

That applying RPS to M{t) is suspect, makes one wonder whether a tighter bound than 

= 1 can be found, if longitudinal causality constraints are properly taken into account. 
Incorporating them would provide a better understanding of holographic result (1.4). In 
Sec. VII we propose another approach — different from the family of RPS models — for the 
computation of entropy in the infinite scattering limit, which again gives = 1. In this 
regard, it is also interesting to note that in the Floquet systems discussed in [19] , there is an 
exact causal light cone for any local operator and u® = 1 is actually achievable, at least for 
a certain class of initial states. But we should note that the systems of [19] do not appear 
to have a continuum limit and thus may not be describable as a relativistic system with a 
local Hamiltonian. 

Finally, let us note that (6.10) is based on tracing over a local Hilbert space at t = 0. The 
tracing is not local at time t. This is a consequence of unitary transformations we performed 
so as to track the trajectories of particles. This is certainly not ideal. In (1 + l)-dimension, 
the equivalent proposal (6.21) motivated from the tsunami picture is based on tracing out a 
local Hilbert space at time t, and thus is conceptually more appealing. But it appears not 
easy to generalize (6.21) to higher dimensions for regions of general shapes as the behavior 
of entanglement tsunamis become complicated at late times. 

From the success of the maximal RPS model in reproducing intricate holographic results 
for multiple intervals in (1 + 1) dimensions, and the discussion of Sec. VIF, it is tempting to 
speculate that a model which incorporates RPS, SSA, and full causality constraints (includ¬ 
ing the longitudinal causality constraints) may go a long way toward describing entanglement 
propagation in interacting systems, and in particular should provide a microscopic physical 
understanding of holographic results. Such a model, however, appears hard to come by from 
the scattering picture. For example, it is not clear how to formulate a precise set of causality 
conditions on the partition of Hilbert space (on which RPS is based). In the next section we 
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use intuition derived from tensor networks to construct an entropy function which coincides 
with holographic results for any number of intervals. 


VII. A NEW MODEL INSPIRED BY TENSOR NETWORKS 

A. Tensor network interpretation of the scattering picture 

We can view the wave function |t^(t)} prepared by the interacting model of Sec. VIB as 
a tensor network, and abandon the idea of associating entropy to the particles ending up in 
the region A. Instead we can regard the time evolution as a quantum circuit of depth t that 
prepares an entangled state from a product state through the action of the unitary scattering 
matrices U introduced in Sec. VI A. It will be easier for us to convert [/ into a four-indexed 
tensor Vj.fcn and to be agnostic about the dimension of the Hilbert space associated to the 
particles scattered, and simply denote it by y; the indices run over i, j, • • • G (l,x). We 
explain our notation in Fig. 14. 

The unitary time evolution prepares us a state 

X 

(t) \i 1 i 2 i 3 ■ ■ ■ ijv) , (7.1) 

{h=i} 

where Ti^i 2 ...iN{^) ‘^an be obtained from contracting Vj.fc/ according to the pattern described 
in Fig. 14. This description provides a convenient interpolation between the free streaming 
and the infinite scattering pictures. In the free streaming case Vi.w = SuSjk^^, while for 
strong scattering we expect V to be random. By tuning V, we should be able to learn how 
the behavior of entanglement spread interpolates between the two. 

The network of Fig. 14 resembles that for global quench described in [2], but there is 
a fundamental difference. Here the vertical direction is physical time, i.e. Fig. 14 is a 
quantum circuit, while in [2] the vertical direction is an auxiliary RG time. Accordingly, 
here the state of the system at a given time t is described by a single slice of the network, 
while there the state is described by the whole network. Nevertheless, in the linear growth 
regime there appears some isomorphism between the network here and that of [2] . It would 
be interesting to understand this better. 

Of course, in this case one has to supply a nontrivial locally entangled initial state. 
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FIG. 14. Top: Explanation of the notation Vij^kh ^tnd an example of the wave function created 
by three particles undergoing two scattering events. The two scattering events can be regarded as 
two unitary gates in a quantum circuit, and the wave function it prepares from the product state 
1000) is represented by the tensor network. The tensor can be obtained by performing one 

contraction on the internal index j. Bottom: The quantum circuit of depth t composed of the 
scattering matrices V prepares \ip{t)) from the product state n^=i|0)Q'- l'^(^)) can be decomposed 
according to the basis \i 1 i 2 i 3 • • •'^w) (7.1), and the coefficient tensor 7iP2^3...^Ar(^) given by the 
tensor network on the figure. 

For tensor networks there exists a bound on the entanglement entropy of a region A\ 

5'^<4utlogx, (7.2) 

where is the length of the minimal cut through the tensor network that separates A 
from A. An example of a minimal cut for two intervals in (1 + 1) dimensions is presented in 
Fig. 15. For large enough % and for a generic unitary scattering matrix V the bound (7.2) is 
expected to be saturated. We then point out that from Fig. 15 it follows immediately that 
in (1 + 1) dimensions our tensor network exactly reproduces the holographic result (6.14) 
for arbitrary number of intervals. In the next section we discuss a higher dimensional model 
inspired by the minimal cut bound in tensor networks. 

We emphasize that the tensor network of Fig. 14 is a purely field theory construct, which 
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FIG. 15. We draw the minimal length cut at time t corresponding to two intervals in (1 + 1) 
dimensions at some intermediate time t. The that cut is not unique due to the discreteness of the 
tensor network. The number of links the cut intersects, f'cut is what appears in (7.2). Note that 
the horizontal section of the cut does not intersect with any links, and it readily follows that £cut 
is proportional to the holographic result (6.14). 

we obtained from the picture of scattering quasiparticles, and a priori it has nothing to 
do with holography. The minimal cut in Fig. 15 is also a tensor network concept. There 
are tantalizing connections with holography though. With physical time replaced by RG 
time, the network here resembles that of [2] which in turns looks like a “nice slice” inside 
an eternal black hole. The minimal cut prescription is also reminiscent of the holographic 
extremal surfaces [20]. 



B. A geometric model inspired by tensor network 

Working with a discrete tensor network in higher dimensions is inconvenient due to the 
breaking of rotational symmetry. Using the physical insight from the tensor network picture, 
we now propose a simple continuum model for the time dependence of the entropy of an 
arbitrary region that satisfies all the physical criteria we are aware of that such a function 
should satisfy. We emphasize that we do not know of any local Hamiltonian that would 
produce the entropy given by this model, and it would be very interesting to find further 
criteria that the entropy function should satisfy that would constrain or rule out the model 
we propose in this section. 

The entropy function after a global quench should satisfy the following geometric require¬ 
ments, some of which we have discussed in previous sections, others were implicit: 
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(a) 

(b) 


(c) 

(d) 

(e) 


(f) 


(g) 


(h) 


The entropy *S'(A) is finite and positive for any spacelike region ^ for t > 0. 

S{A) is the same for any Cauchy surface for A, or phrased in another way, *S'(74) is 
a function of the domain of dependence of A. This is a requirement for all Lorentz 
invariant theories. 

Take any Cauchy slice of the spacetime that contains A, and define A as the comple¬ 
ment of A on this slice. Then 5'(74) = >S'(74). 

S{A) is invariant under translations and rotations of A. This condition follows from 
homogeneity and isotropy of the state. 

For any space or lightlike region A the scaling relation S{A^l’^^) = X'^~^S{A) holds, 
where we defined the (in time and space) scaled region = {\x\x E A} and 

A > 0. This relation is just the generalization of (4.6) to regions which do not lie on 
a constant time slice. 


For any region A lying on a constant t slice define the spatially scaled region at the 
same time by A^^^ = {(t,Xx) \ (t,x) E A}. We have the volume law for small regions, 


lim 

A^-O 


voi(A$;"^) 


= s. 


(7.3) 


For any region A lying on a constant t slice define the time scaled region by ^4^^ = 
{(At, x) I (t, x) E A}. We have the tsunami law for small times, 


lim 

A->0 


area(9A^dj ;\ 


SVe . 


(7.4) 


Strong subadditivity S{A) + S{B) > S{An B) + S{AU B) applies for any two regions 
A and B on the same Cauchy surface. 


Using some insight from the tensor network picture, below we will produce a family 
of entropy functions *S'(74) — including one with ve = I — with these properties in any 
dimension. Of course, the free streaming models of Sec. IV also obey all these properties, 
but have ve < 
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For the construction let us introduce the set of all [d— l)-dimensional surfaces on 
> 0, that are smooth almost everywhere^® and have a normal vector n = (n°,n), which 
satisfies in all points where the normal is well defined 


'n^ 


< 


1 


\n 


01 — 


0 < a < oo . 


a 


(7.5) 


That is, the slope of the tangent vectors of the surface is bounded above by a; a = 1 means 
that the maximal slope is given by the speed of light. 

For any S G $0 let us define V (S) as the volume of S using the metric ds^ = dx'^, which 
is degenerate. That is, V(S) is the volume of the spatial projection of E. Then, let us 
consider the following entropy function for an arbitrary spatial region 


Sa{A) = s min y(E). (7.6) 

For this formula we can consider surfaces S that hit the t = 0 boundary of the spacetime 
and end there. Equivalently, we can think that once they hit t = 0 on the boundary, they 
run on the t = 0 plane, but this part of the surface does not have any volume. Because for 
only the spatial volume counts, a large number of surfaces give the same entropy. The 
minimal surfaces E can be thought of as the analog of a minimal cut in the tensor network; 
as explained in the caption of Fig. 15, the minimal cut is not unique either. 

It is immediate that items (a), (b), (c), (d), (e), and (f) above are satisfied for any a < 00. 
Regarding the tsunami velocity for small times, we see that the least volume will be given 
by a surface that runs to the past of A as fast as it can subject to the constraints on the 
slope to end on the t = 0 boundary. This gives 


VE = a. (7.7) 

From the arguments of Sec. II it is clear that ve > 1 violates SSA, which is condition (h) 
above. This then restricts the range of a to 0 < a < 1, and we will show below that for any 
a in this range (h) is satisfied. In particular, for a = 1 we have a model that gives ve = 1 
and satisfies all the conditions we listed above. 

We now prove that SSA holds for 0 < a < 1, but is violated in some circumstances for 
a > 1. Let us take a Cauchy slice containing two overlapping regions A and B, and for 

That is, they can contain singular codimension-l sets, where the tangent is not defined. 


49 




0 < a < 1 use the degeneracy of the minimal surfaces to push S^ and to the past so that 
they necessarily intersect, see Fig. 16. We can reinterpret the two surfaces as corresponding 
to AU B and A H B. They may not be minimal, but minimal surfaces corresponding to 
AU B and An B will just decrease the entropy, so SSA follows. These steps are identical 
to the proof of SSA in holography for static situations [21]. 



FIG. 16. Left: Minimal surfaces drawn by dashed lines corresponding to a strip A. Because ds‘^ 
is a degenerate metric only the spatial projection of the surfaces contribute to their volume F(T). 
Here we chose o < 1, so the surfaces from are allowed to be steeper than the light cone. Right: 
Proof of SSA for the case 0 < o < 1. We choose A to consist of a constant time and a lightlike 
strip (drawn by blue), while R is a union of two lightlike strips (drawn by green). (The proof works 
for arbitrary A and B on the same Cauchy slice.) Example minimal surfaces corresponding to A 
and B are drawn by dashed lines and colored in the same colors as A and B respectively. We can 
reinterpret the intersecting minimal surfaces as surfaces ending on AdB (drawn by red) and AuB 
(not drawn). Note that in the present case the red surface is minimal, while the one corresponding 
to AU B may or may not be minimal depending on the position of the t = 0 boundary. 

The above proof fails for a > 1, because in some cases the minimal surfaces corresponding 
to A and B cannot be made to intersect. SSA is violated in the setup of Fig. 17, as the 
entropies for A, S, and AUB are all proportional to their spatial projection, while S(AnB) 
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is larger than the spatial projection of A H B due to the turning point. Note that for a 
lightlike (or sufficiently boosted) strip the minimal surface is necessarily cuspy in this case, 
any smoothing would bring the surface out of the set we allowed for codimension-1 
singularities on S in our assumptions. 

This model for u® = 1 {a = 1) in (1-1-1) dimensions gives exactly the prescription of 
the tensor network discussed in the last subsection, and it coincides with the holographic 
result. For higher dimensions it again leads to = 1. Note that in this model in higher 
dimensions the minimal surface that gives ve = f corresponds to the situation where entan¬ 
glement essentially only propagates along one direction (i.e. the direction perpendicular to 
the boundary), which is clearly a bit peculiar. This aspect is similar to the behavior of the 
maximal RPS model discussed in the last section. It would be nice to understand whether 
we are missing some more subtle isotropy constraint. Without using the isotropy of the 
system, one cannot hope to prove a more restrictive upper bound than < 1, as many 
decoupled (1 -|- l)-dimensional CFTs placed next to each other satisfy all our assumptions, 
and trivially yield ve = 1. Another question we leave for future research is how the present 
model compares with the holographic results once we set a = 



FIG. 17. Failure of SSA for a > 1. We use the same regions A and B as on Fig. 16. However, 
the minimal surface corresponding to A now cannot be pushed further to the past, while the one 
corresponding to B cannot be pushed further to the future. Thus, minimal surfaces corresponding 
to A and B cannot intersect, and the geometric proof of SSA breaks down. SSA is violated because 
the minimal surfaces corresponding to AnH have larger volume than the spatial projection of AnH. 

Finally, we note that the family of models we constrncted in this section gives for a region 
lying on a constant time slice <S'(A) < svol(A) for all times, and can give rise to nonzero 
mutual information at intermediate times. 
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C. Holographic theories maximize entanglement spread in (1 + 1) dimensions 


The tensor network model seems to include by construction the heuristic idea of a “max¬ 
imal spread of entanglement”. If we think that strong interactions are effectively taken into 
account by random unitaries in the vertices of the tensor network, the model would try to 
produce locally a random pure state as fast as possible given the causality constraints. 

In fact, we can prove that the holographic entropy formula Sh{A) (6.14) gives the absolute 
maximum for the entropy *S'(^), for any number of intervals, and any global quench in 
relativistic theories in (1 -|- 1) dimensions. 

In order to show this, let us take regions at a fixed time t and start with a single interval 
of size r. From SSA we have S"{r) < 0 (for any translationally invariant state). From 
the physical assumptions about the global quench we also have that S{r) ~ sr for small 
r, as small intervals saturate at the volume law, and that the entropy saturates to an r 
independent constant S{r) ~ 2svEt for large r, as large intervals are in the linear regime. It 
is immediate that the holographic Snir) = smin(r, 2vEt) is the maximum concave function 
with these two asymptotic behaviors. Here ve = 1, but we write it explicitly for convenience. 

For n intervals the proof is by induction. We start by assuming that Sh{B) is the 
maximum possible entropy for any set B with less than n intervals. 

For n intervals let us write A = {(ai, bi), (02,62), • • •, (dn, bn)}- Recall that the geodesics 
describing the bulk extremal surface for this region in the holographic setup all lie in some 
spacelike surface, and cannot cross each other [22] (otherwise another surface exists with 
smaller area).^° It is not difficult to realize^^ that there must be a geodesic in the extremal 
surface for A that either joins the two end-points of an interval from A, say Ij — (aj,bj), 
or there is at least one geodesic that joins the two consecutive points bk,ak+i for some 
k. This last possibility is equivalent to joining the two end points in an interval from A, 
which itself is a union of intervals and two half lines. As these two cases are completely 
analogous, without loss of generality let us restrict our attention to the first case, i.e., that 
Ij = {aj,bj) determines a geodesic in the extremal surface. Then the rest of the geodesics 

That Sh{A) gives the absolute maximum for the entropy can be proven directly from the formula (6.14) 

without referring to geodesics, but perhaps it is easier to visualize the proof presented here. 

This can also be proved by induction: a given geodesic splits the set of geodesics in two, to the ones 

outside and inside of it, and consequently the problem is mapped to the same problem with less intervals. 
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give the extremal surface iov B = A — Ij = A H Ij. We have 

Sh{A) = Snilj) + Sh{B) > S{Ij) + SiB ), (7.8) 

where in this last inequality we used the induction hypothesis. Now, by subadditivity 

S(I,) + S(B) > S(A) , (7.9) 

and this gives the proposed inequality 

SHiA)>S{A). (7.10) 

We note that the assumptions about the function S{A) are minimal: we only used SSA for 
regions on constant time slices, the volume law for small regions and saturation of entropy 
for large regions. In this sense the proof can be carried over to nonrelativistic situations, 
or to systems with < 1. The proof can easily be generalized to regions not lying on a 
constant time slice in the relativistic case. 

From another perspective, this result also sheds different light on holographic entangle¬ 
ment entropy itself. For example, essentially the same proof shows that the holographic 
entanglement entropy given by minimal geodesics in pure AdSs space is an absolute maxi¬ 
mum for the vacuum entanglement entropy of any region and any (1 -|- l)-dimensional CFT 
with the same central charge. 


VIII. TSUNAMI VELOCITY AND MUTUAL INFORMATION 


In this section we further examine the relation (2.9), which we copy here for convenience 
(see also Fig. 2) 




( 8 . 1 ) 


in the context of free streaming and interaction RPS models discussed in earlier sections, 
which sheds new insight onto the values of 

Let us first consider (1 -|- 1) dimensions. In a free propagation system with all particles 
traveling at the speed of light, from plot (a) of Fig. 18, it is clear that there is no entanglement 
between W{t) and X (see caption). This is consistent with u® = 1 in such a model. There 


22 


We need the condition of equal central charge to have the same entropy for a single interval. 
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is, however, entanglement between W{t) and X even in a free propagation system if there 
are also particles which can travel smaller than the speed of light. See plot (b) of Fig. 18. 
This means that whenever there is propagation of entanglement inside the light cone we 
expect ve < i- 

In a system with only particles traveling at the speed of light, interactions will generically 
generate entanglement between W{t) and X, as indicated in plot (c) of Fig. 18. Thus in 
(1 + 1) dimensions, interactions slow down ve compared with the free propagation at the 
speed of light, and for a general interacting theory we should have ve < 1. This conclusion 
is consistent with ve = 1 from CFT [1], as well as holographic calculations [16], where the 
results apply in the “infinite scattering limit” R,t ^ In other words, in these theories, 
we expect 

MsW = 1 - ^/(*/-R) + ■ ■ ■ (8-2) 

with / a positive function. For holographic theories, the function / can be read from (5.38) 
of the second reference of [3] and is indeed positive. [17] also found that interactions that 
break the scaling symmetry slow down the tsunami. It can also be readily checked that in 
the maximal RPS model I{W{t),X) = 0. 

The story in higher dimensions is rather different. Even in a free propagation system 
with all particles traveling at the speed of light, there is entanglement between W{t) and X, 
which can be immediately seen as follows. Particles from a point x with nonzero velocity 
in directions perpendicular to xi will have a velocity smaller than the speed of light in the 
Xi direction. In other words, when projected to the Xi direction, these particles propagate 
“inside” the light cone. Thus from plot (b) of Fig. 18 they will generate nonzero I{W(t),X). 
This is perfectly consistent with our earlier discussion that u® < 1 for higher dimensional 
free propagation models. In fact by computing directly X) and then using (2.9) 

gives an independent derivation of ve, which by consistency should agree with our earlier 
expression (4.27). Below we will check this is indeed the case. 

With already nontrivial entanglement between W(t) and X in free propagations, interac¬ 
tions can increase u® if they can reduce entanglement between W (t) and X. It appears hard 
to imagine that the entanglement can be reduced completely to zero. But at the moment 
we do not have any definite argument to exclude the possibility that in an isotropic system 
there exists some limit, in which = 1 can be approached like in (8.2). 

Finally, as promised earlier we show that ve derived from (2.9) agrees with (4.27) for a 
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FIG. 18. Entanglement of W{t) and X in (1 + 1) dimensions. The definition of W{t) and X are 
the same as those in Fig. 2. Plot (a): In a free propagation model with particles traveling at speed 
of light, only particles from region labelled by can contribute to S{X). Clearly there is no 
entanglement between W{t) and X, as no particles from the same point can both reach X and 
W{t). Plot (b): In a free propagation model which contains also particles traveling smaller than 
speed of light (i.e. inside the light cone), there is entanglement between W{t) and X. The region 
M{t) which can contribute to S{X) is also much larger. Plot (c): Even with particles only traveling 
at the speed of light, scatterings can also generate entanglement between W{t) and X. Suppose 
initially 1 and 2 are entangled. Scattering between 2 and 3 will generate entanglement between 1 
and 3, thus leading to entanglement between W{t) and X. The region which can contribute 
to S{X) is also much larger than that of (a) and is the same as that of (b). 
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general free propagating model. Consider a light cone from a point of distance y from the 
entangling surface (see e.g. Fig. 19), as in Fig. 6. Positive y is outside E, negative is inside 
E. A light cone centered at y intersects the regions W{t) and W(t)\J X in spherical caps 
with opening angles 

9 = arccos , (8.3) 

9 + 56 = arccos ( -—^ ) , (8.4) 

\t-\- 5t J 

respectively. The corresponding spherical caps have normalized area ^ and ^+(5^ respectively. 



FIG. 19. W{t)UX is drawn by blue, the light cone by purple, and the entangling surface by green. 
The left figure is for y > 0, while the right figure depicts y < 0. 


In complete analogy with the entropy, the mutual information can be calculated for every 
light cone independently, as in (4.1). It is given by a formula analogous to (4.28) 

I(W{t),X) = Aj: j dy {fice.p[^iy/t)] +SeqS^iy/t) - + S^{y/t)]), (8.5) 

where we used that the small annulus region that we get from the difference of the spherical 
caps has to obey (4.4). On Fig. 19 this region is the two small arcs between the blue and 
green planes. We series expand (8.5) to get 

I{W{t),X) = dx <^^(^), (8.6) 

where we introduced x = y/t. We calculate the first term using the expressions: 


5^{x) 

59{x) 


^ 9{x) 59{x), 

^d-2 

5t 1 — cos 9(x) 
t sin0(a:) 


(87) 

( 8 . 8 ) 
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(8.9) 


Using these expressions we obtain after the change of variables dx = dcos9: 


Seq dx 5^{x) = St . 


-1 


For the second term in (8.6) we use the chain rule 



( 8 . 10 ) 


to write: 



( 8 . 11 ) 


( 8 . 12 ) 

(8.13) 


= -Seq St Ve 1 


where in the second line we did a partial integration and used /rcap(O) = /Ucap(l) = 0, while 
in the third line we used (4.27). Combining the two terms then gives 


I{W{t),X) = s,^Aj:5t{l-VE) 


(8.14) 


which reproduces the identity (2.9). 
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Appendix A: Concavity of the spherical cap entropy function 

Eirst, we consider a situation in d = 3, where we apply (4.3) to two regions A, B shown 
in Fig. 20, with all regions A, B, A Cl B, A U B singly connected and included in half of the 
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light cone. We then hnd 


/^(Ca) + /^(Cb) > l^iUnB) + KCaub) 

which immediately leads to 


(Al) 


< 0, (A2) 

if we take = C; Cb = C) Caob = C ~ e, ^nuB = C + ^ with e —)■ 0. That is, fi{^) is a concave 
function. 


A B 



FIG. 20. Regions A, R, A n R, A n R on the light cone drawn by red. 


Second, we consider spherical caps in d > 3, where the geometry is more complicated. 
The simple proof of the concavity of Pcap in d = 3 is similar to the entropic proof [23] of 
the C-theorem in d = 2, whereas the proof in d > 3 below uses the more complex method 
of [24] that was originally developed to prove the F-theorem in d = 3. 

We are going to use strong subadditivity to prove the concavity of the function /icap(0 
as a function of the normalized area Strong subadditivity (4.3) for two spherical caps 
A and B involves the intersection and union of these regions, which are not spherical caps 
anymore. In order to have an inequality containing only spherical caps we can follow the 
procedure used in [24]. We use a large number N of spherical caps X*, i = 1,..., X, which 
are copies of one spherical cap rotated around a point in the sphere, see hgure 21. Strong 
subadditivity leads to 

N 

^ S{UiXi) + n Xj)) + S{U{ijk}{Xi nx,- nx,)) +... + R(n,w) • (A3) 

i=l 

The regions on the right hand side of this inequality are not spherical caps, but as shown 
in hgure 21 they approach to spherical caps in the limit of large In this limit (A3) 

We know the wiggly caps of figure 21 have an entropy which converges to the one of the spherical cap 
because, as we have shown in Sec. V, variations of entropy are bounded by Seq times variations of area. 


58 







FIG. 21. N spherical caps which are copies of a single spherical cap rotated around a point different 
from its center (shown in blue). The unit sphere is shown in red and the region with black contour 
is one of the regions appearing in the right hand side of (A3). 

will be converted into an inequality involving an integral of spherical caps with sizes varying 
between a maximum ^max and a minimum ^min (see figure 21). Following [24] we can then 
take — ^min = ^ and expand the inequality for small e to get an inequality for /icap(0 
and its derivatives for a single 

Though we can follow these same steps here, we can obtain directly the final result 
without doing the explicit calculations by arguing as follows. As a result of this procedure 
we should get a differential inequality for /Ucap(0- Because strong subadditivity is linear in 
the entropy and involves four regions, it will always lead to linear differential inequalities 
containing at most second derivatives of the entropies. Hence we should get 

+ /2(e/‘c.p(0 + /3K)fcpK) < 0. (A4) 

Now, it is evident that the constant function S{A) = const, is a solution of the strong 
subadditive equation S{A) + S{B) = S{An B) + S{A[J B) rather than inequality. Hence, 
a constant /icap(0 = const, must be a solution of (A4) with equality rather than inequality. 
This is only possible if /3(^) = 0. In the same way, the area function 5'(A) = is a solution 
of strong subadditivity with equality for any regions. Hence /rcap(0 = C should be a solution 
of (A4) with equality to zero. This implies /2(C) = 0. Therefore we get that /r'ap(0 always 
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either positive or negative for ^ according to the sign of /i(^). But we know the examples 
discussed in this paper all have /icap(0 — 0- This then implies /i(^) is positive, and in turn 
the general inequality 

(•IpK) < 0 (AS) 

for any entropy function on the sphere. 

Appendix B: Selected consequences of ballistic propagation of entanglement 

In this appendix we study in more detail various aspects of entanglement propagation 
using the three examples discussed in Sec. IV B. We will mostly use spherical regions for 
illustration, although in some cases we make general remarks valid for all measures and 
general shapes. More specifically, we will investigate the following issues: 

• full time evolution and the entanglement rate (1.5) for simple shapes such as the 
sphere, 

• for the EPR example, one can further define an entanglement density which can be 
used to better visualize the spread of entanglement, 

• saturation time for generic shapes, 

• finite volume effects, 

• mutual information. 

Along the way we also compare the results obtained here for ballistic propagation with 
those of holographic systems. 

1. Time evolution for simple shapes 

a. The d — 2 interval 

Let us first briefly review the results in d = 2 [1]. The entangling region reduces to an 
interval and the geometry of the problem is presented in fig. 22. On the figure the green 

It is important though that we are assuming a finite and smooth entropy function (4.4), in contrast to 
the vacuum contribution to the entropy, which would give a divergent area law piece. 
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region is Afit), defined as the region of the space where he the centers of the fight cones 
that have non-empty intersection with E at time t. In one spatial dimension, because the 
intersection of a light cone with the region E is a point, all entanglement measures are 
equivalent. From the perspective of quasiparticle propagation, since there are only two 
directions and all particles in a direction from a single point propagate side by side, no 
matter how entanglement is distributed, the entanglement spread reduces effectively to that 
of the EPR example. 

We take the width of the interval to be 2R. From Fig. 22 we conclude that the entangle¬ 
ment entropy has the time dependence: 


s jit {t < R) 

^ [dR {t>R), 


(Bl) 


where s is given by (4.9). We divided s by 2 as n/2 is the density of quasi-particle pairs. 
Both the slope of the linear growth and the saturation time agrees with those found from 
direct field theory calculation [1] and those obtained from holography [16]. 


2R 



2R 



2R 



FIG. 22. d — 2 spacetime diagram. The first figure is for t < R, the second for the saturation 
time t — tg — R, and the third after t > R. The length of the green intervals determines the 
entanglement entropy. 
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b. Full time evolution for a sphere 


Let us now take spherical entangling surfaces of radius R, for which the high symmetry 
enables an analytic treatment of the full time evolution. For a spherical region, the intersec¬ 
tion A of any light cone with the region is always a spherical cap, and the EPR, and RPS 
measures all contribute to ii[A] = p(min(^^, ^^)). We can conveniently treat all measures 
at once by working with: 




(B2) 


where m = 1 gives the result for the EPR and RPS measures, while m > 1 corresponds to 
GHZ blocks. The absolute value comes from combining the relation = 1 with (4.13) 
valid when < 1/2. 

To apply (4.1) we need to work out first the region N'sd -2 (t), which is explained in Pig. 23. 
We then find that 


p pt-\-R 

S(t) = / d‘^~^r iJ,capi^ir)) = oJd -2 / dr iJ,capi^ir)), 

djVsd-2it) J\t-R\ 


(B3) 


with 


min(^(r), l-C(r)) 


2 'V 2 ’2/’ 


{R + r + t){R + r — t){R — r + t){R — r — t) 
4^2^2 


(B4) 

where we used the formula for the area of a spherical cap, and Iz{a,b) is the regularized 
incomplete beta function. We did not find a way to perform this integral for general d, m, 
so we present some examples. 

Por the EPR case m = 1, we find (with r = t/R) 


Sd=z{r)/R^ 


Sd=A{T)/R^ 


Sd=^{r)/R^ 


{ 2s [rVT — + arcsinr] (r < 1) 

ITS ('T > 1), 


127^5 (r - y) (^ < 1) 
(r>l), 

J arcsin 


(r < 1) 

(r > 1) . 


(B5) 


(B6) 


(B7) 
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FIG. 23. The region Msd- 2 {t) for d — ?>. The left two figures are for t < i? and the right ones are 
for t > R. The top figure in each column gives two “critical” light cones (dashed circles) which 
just touch the entangling surface (solid circle). They can be used to determine the boundary of 
Nsd- 2 {t)^ which is the shaded green region in the bottom row. For t — R we would get a filled 
green disk of radius 2R. Note that if we restrict to the x-axis we get back the d = 2 picture of 
fig. 22. 

Numerical plots for various m in d = 3,4 for the GHZ measure are given in Fig. 24. Note 
that while for EPR and random state measure, the saturation time is always G = i?, for 
GHZ the saturation time is infinite. It can also be readily checked that for early times 

S (t) = sbJd-2 t + ... (B8) 

and at late time the entropy Ssd- 2 {t —?• oo) = s consistent with our general 

discussion in Sec. IV C 1. 

c. Strip in d = 3 

We consider a strip of width o and length L in d = 3 as another simple example. Our 
discussion will be less detailed, than for the sphere case. For early times, t < a/2 the time 
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S(t)/(s^R^) S(t)/(s^R^) 



FIG. 24. Top row: Full time evolution for a sphere in d = 3,4 for m = 1,..., 5, with larger m 
corresponding to darker color. Bottom row: Rate of entanglement growth Digd -2 (1.5) obtained 
by taking the time derivative of the functions in the top row. Note that '^gd -2 is monotonically 
decreasing for all m, and hence is bounded by the tsunami velocity 


evolution is exactly linear, while for t > a/2 the width of the strip starts to play a role. 
The time dependence becomes more complicated as for some of the light cones Lstrip(x;t) 
becomes a disconnected region, and we have to use (4.14). The result for early time valid 
for arbitrary measure is 


Sanp(T)/(aL) -2 sveT (t < 1 / 2 ), 


(B9) 


with T = t/a. For the EPR model, we get for later times 


2s 


'S'strip (^)/(®.b) 


TT 


2r — \/ 4t 2 — 1 + arccos — 

2r 


(t > 1 / 2 ) 


(BIO) 


For the case of the strip the RPS model is not equivalent to the EPR model. The results 
for the RPS, the EPR, and GHZ block models are plotted in Fig. 25. It is remarkable that 
the measures originating from a quasiparticle picture give an inhnite saturation time, while 
the RPS model saturates in r = 1. 
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S(T)/(;yeqLa) 



FIG. 25. Full time evolution for a strip in d — 3 for the GHZ block model with m = 1,..., 5, with 
larger m corresponding to darker shade of orange. The RPS model gives a different time evolntion 
than the EPR model and is plotted with blue. 


d. Two disks and strips in d = 3 

We now briefly examine the entanglement entropy of two separated disks and strips. We 
use the RPS and EPR measures. These examples are complicated enough geometrically to 
give interesting features in the time evolution: we show on Fig. 26 that the entanglement 
growth rate (1.5) turns negative for them at some intermediate time for the EPR measure. 

This phenomenon is the easiest to understand in the case of the two disks. We chose their 
separation so that at some intermediate time the entropy saturates. However, at some later 
time some of the quasiparticle pairs that originated from between the two disks do not give 
entanglement with the outside, as one can ends up in one disk, the other in the other disk. 
Because the system has already reached the saturation value for the entropy at intermediate 
times, the entropy must go down. This resonant effect only last for a finite time, and finally 
we get complete saturation. We discuss aspects of saturation in Appendix B3. Similar 
discussions can be found in [1, 8-10]. 

We included the example of two strips, as it also displays negative entanglement growth 
rate, and this geometry will also play a role in our discussion of mutual information in 
Appendix B 5. In Fig. 26 we show a setup, where both the EPR and RPS measures exhibit 
a dip in the entanglement entropy. 
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FIG. 26. Full time evolution for two disks and strips in d = 3. The confignrations are shown in 
the top row. The radins of the disks is i?, and they are separated, by a distance 2R. The width 
of the strips is a, and they are separated by a distance I — a. r is the dimensionless time equal to 
t/R and t/a respectively. For the case of two disks we have drawn an EPR pair that contributes 
to the resonant effect shown in the bottom row. In the bottom row the time evolution is plotted 
for the RPS (blue) and the EPR (orange) models. 


2. Entanglement density 

For the case of EPR pairs we can introduce a local, and thus more refined, measure of 
entanglement: entanglement density. The entanglement density p^(x,t) at a given point 
X E A inside S is defined as u times the density of quasiparticles whose entangled partners 
he outside S. Recall that u was introduced below (4.9). It then immediately follows that 

ST.it) = I d^~^xpT{xA) ■ (Bll) 

Ja 

PTix, t) can be readily worked out as follows. In the EPR example, two point x, y are are 
only entangled for one moment, when their distance is exactly 2t. One can then introduce 
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an entanglement “correlation function' 


s{x,y) = 


5{\x-y\ -2t) . 


UJd-2 

The normalization of the above function is determined by requiring 


(B12) 


IT ^ys(x, y) = s = nu, 


(B13) 


which is the total amount of entanglement x has with the full space. The entanglement 
density (T, t) can then be obtained by integrating the above expression over all y that lies 
in A, i.e. the region outside E, 


P2 (f,t)= [ ^ys{x,y) = - [ dy 5{\x-y\-2t) . (B14) 

JA ( 2^1 ^d-2 Ja 

The above expression can be easily described in words: draw a sphere of radius 2t around 
X, and calculate the portion of the surface that falls outside E. The simple intuitive picture 
this is that we are counting the quasiparticles at x that have their partners outside, which 
all lie on the sphere of radius 2t drawn around the point. 

For E a sphere, the integral in (B14) can be readily performed. Actually, we have already 
performed this calculation in Appendix Bib, so all we have to do is to replace t ^ 2t 
in (B4). For d = 3 we get 


J R^-r‘^-At^ \ 
Art ) 


p{r) = s < 


0 


(|i? — 2t\<r<R) 
{Q <r < R — 2t) 


1 {0 <r <2t - R), 


(B15) 


where the last two cases can only happen before and after t = R/2. The entanglement 
entropy saturates when p(r) reaches s everywhere inside the sphere. It is curious that at 
the center of the sphere r = 0, the density jumps from p{r = 0) = 0 to saturation value s at 
t = R/2. We plot (B15) for various times in Fig. 27. 

The entanglement density can be used to give a precise definition of the entanglement 
tsunami introduced in [3] : we define the tsunami wave front as the boundary between regions 
of p = 0 and p 7^ 0. From (B14), one can immediately conclude that under such a definition, 
locally the wave front should progress at the speed of 2, i.e. twice the speed of light. For E 
being a sphere, the wave front can be visualized from the lower half of the plots in Fig. 27, 
with the wave front reaching r = 0 at t = R/2. Note that when the wave comes in, the region 
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covered by the wave is only partially entangled with outside, i.e. with a value smaller than 
the equilibrium value s. The curves in the upper half in Fig. 27 suggest there is a “reflected 
wave”, whose wave front can be defined as the boundary between the region which has 
reached the equilibrium value s and the region which has not. This reflected wave starts 
Ski t = R/2 from the center and moves at speed 2 outwards, reaching S at t = i?. In this 
example, the linear growth (4.26) and the associated ve can be considered as an average 
effect. The picture here is very different from that proposed in [3] for strongly coupled 
systems, where the region covered by the tsunami wave will already have reached their 
equilibrium value. The difference may be due to that in free theory as we are considering 
here, there is actually no dynamical process of equilibration. We should also keep in mind, 
as we will also elaborate below, that generically there is no unambiguous definition for an 
entanglement density. 


P/^eq 



FIG. 27. Entanglement density as a function of radial distance. Different lines show to the density 
profiles at different times, and later times correspond to darker colors. The straight line in the 
middle is for t = R/2. The curves below t — R/2 are for t < R/2 with time increasing from right 
to left. The curves above t — R/2 are for t > R/2 with time increasing from left to right. 


The entanglement density discussed above is specific to the EPR example. When there 
is multipartite entanglement, such a definition does not appear to exist even assuming a 
quasiparticle picture. The basic reason is that with multipartite entanglement, one could 
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no longer localize entanglement to a point. This can be readily seen from the GHZ example 
illustrated in Fig. 28. We discuss how interactions change the perspective on the tsunami 
wave front in Sec. VIC. 



FIG. 28. Explanation of why the entanglement density cannot be defined for GHZ blocks. In the 
figure we take m = 2. For the left plot we could localize the entanglement to the orange point 
inside S, but in situations depicted in the middle and right plots, the entanglement among the 
block can no longer be localized to the orange point. 


3. Saturation time 

In Sec. IV C we showed that the entanglement entropy has an equilibrium value S'eq = sV^. 
In Appendix B 1 b we saw that for a sphere of radius R, the EPR model has a finite saturation 
time given by G = R, while for GHZ block example the saturation time is infinite. In this 
subsection we make some general remarks on the saturation time for general shapes. 

Let us first consider the EPR example. From (B14) we can immediately conclude that 
for any S the saturation time tg equals half of the largest distance between two points on S. 
If t > tg the entanglement density at any point inside S is p{x) = s, because the Dirac delta 
in (B14) has support on a circle of radius 2t around x, which now completely lies outside 
E; see Fig. 29. Phrased slightly differently, at time t, all quasiparticles which entangle with 
those at x he on the circle of radius 2t centered at x, and the total number of them is n, 
the particle density at x (which is in fact the same everywhere due to homogeneity). When 
this circle lies completely outside S, all of them contribute to S's and the entropy does not 
change with time. The saturation time is thus tg = R for spherical S and tg = oo for any 
non-compact region. In particular, tg is infinite for a strip. 

The RPS measure gives a rather different saturation behavior from the EPR model as 
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FIG. 29. The figure on the left is for t < tg, the right one for t > tg. The purple circles are circles 
of radius 2t centered at some point in E. At time t, all quasiparticles which entangle with those at 
a point X lie on such a circle centered at x. For t > tg all such circles completely he outside of E. 

evidenced by the finite saturation time even for a non-compact shape; for the example 
of the strip see Fig. 25. The criterion for saturation in the RPS model is that after the 
saturation time tg there should not exist any light cone, whose intersection with region A 
has a normalized area grater than 1/2. 

The results of the EPR and RPS measures should be contrasted with holographic systems, 
for which tg is finite for a strip and given by 

R 

,,holo ’ 

where 2R is the width of the strip and u® is the “tsunami velocity” which appears in (1.4). 
For holographic systems, the saturation time for a sphere is given by 

R d 

tg = — > R, ce = \ vT-i—^ <1 for d>2, (B17) 

ce y 2 {d- 1) 

where we have quoted the value of ce for a neutral system. We thus see that for a strip, 
holographic systems saturate faster than the EPR model (tg = oo) and the RPS model 
{tg = 2R, as can be calculated from the above criterion or read from Eig. 25), while for a 
sphere holographic systems saturate slower than the free streaming models that have tg = R. 
Note that the velocity ce has also appeared in [25] as the “expansion” velocity of the the 
time evolution of a local operator in a thermal state. 
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For a general measure from our discussion in Sec. IV C 1 we expect that in the t ^ oo 
limit, there are generically 1/t corrections to the leading behavior (4.22). This implies 
that the saturation time is generically infinite. As an example, let us consider S a sphere, 
for which case the shell in Fig. 5 is also spherical. One can then compute the subleading 
corrections using (B3) by expanding /rcap(0 higher orders in 

/^cap(0 = SeqC + 02^^ + + • • • , (B18) 

where we used (4.23). We will not go into details here, except to mention that whenever 
the nonlinear term in (B18) are non-vanishing one gets subleading corrections in 1/t of the 
form, 

Sit) = ■ (B19) 

n 

Thus all measures for which there are nonlinear terms in (B18) have infinite saturation time. 

4. Finite volume effects 

So far our discussion assumed the system has an infinite volume. If the system has a 
finite volume, the large time behavior of the entanglement entropy will be modified when 
carriers of entanglement can explore the whole volume. 


o. EPR pairs and GHZ blocks 


Let us first consider the EPR and GHZ examples which can be treated in a unified manner. 
For a system with a finite volume, after a long time there will be no correlation between the 
positions of the particles that originated from the same point. (This assumptions should 
hold true except for resonant situations in special geometries.) Then we have a constant 
density n of quasiparticles, and a total of nVystem of them, randomly distributed. 

We get entanglement except in cases, when quasiparticles originating from one point are 
all inside or outside S, hence 


S'Ejit ^ Oo) 2^System 


Rs A 

System / 


2m 






system 


2m 


(B20) 


This expression is manifestly symmetric under Vs —t Fsystem — Ffi, hence the requirement 
that the entropy of a region and its complement is equal in a pure state is satisfied. Another 
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consistency check is that in the infinite i4ystem limit we get 5 'e = sVs- The maximum 
of (B20) is achieved for Vi]/f4ystem = 1/2, 




1 - 


1 


22m—1 


(B21) 


In the EPR (m = 1) case the resulting expression is 


Sj:{t oo) = 


14 


14 (^system - 14 ) . 


(B22) 


system 


The maximum entanglement is 


_ ST/ 

max 4'^system- 


b. Random pure state measure 

According to the discussion in Sec. Ill random pure states are expected to give: 

ST.(t^Oo)=S min (14, 14ystem “ 14) • (B23) 

This result is consistent with (4.16); if we followed the time evolution of fight cones in a 
compact geometry, a fight cone would become a curve densely filling the whole volume 
14ystem- The maximum entropy is reached again at 14/14ystem = 1/2, and its value is 
‘^^Imax ~ flsystem, twicc the valuc of the EPR pair model. Equation (B23) is of the form 
expected from a holographic system. 


5. Mutual Information 

We now consider the qualitative behavior of the time dependence of mutual information 
I (A, B) between two regions A and B. 


a. EPR pairs 

In the EPR model a pair contributes to the mutual information I{A, B) = S{A) + S{B) — 
S{AB) for regions A and B, if one member of the pair is inside A and the other is in B. Let 
the smallest distance between the two regions be Lmin, and the largest Lmax- (If one of the 
regions is non-compact L^ax is infinite.) Because of back to back propagation, the mutual 
information is nonzero only for times 

^ ^ ^ _ (B24) 

2 2 
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If the whole system is compact, the mutual information will tend asymptotically to 

I{A,B) = -^VaVb, (B25) 

^ system 

where we have used (B22). 

An interesting feature of the EPR example is that the corresponding mutual information 
is extensive, i.e. 

= or Jepr(A,S) + /epr(A,C') = /epr(A,BC'), (B26) 

where the tripartite information is defined by 

h{A, B, C) = I{A, B) + /(A, C) - I{A, BC) 

= S{A) + SiB) + S{C) - S{AB) - S{BC) - S{AC) + S{ABC) . (B27) 

To see this let us consider the contribution of an EPR pair to Is (A, B, C) for some regions 
A, B, C. Without loss of generality let us assume that neither of the quasiparticles is in C. 
This implies 

S{C) = 0, S{AC) = S{A) S{BC) = S{B) S{ABC) = S{AB) , (B28) 

from which (B26) immediately follows. 

b. GHZ blocks 

For GHZ, /(A, B) becomes non-zero at the time as in the EPR case. Here, however, 
for any time we can always find a light cone that intersects both A and B, and we can put 
some quasiparticles on these intersections without violating momentum conservation. This 
leads to a non-vanishing I[A, B) for all times. As t —>■ oo, as the normalized volume of 
intersection of a light cone with A and B will necessarily go to zero, therefore /(A, B) will 
asymptote to zero. This phenomenon is reminiscent of the differences in saturation between 
the GHZ and EPR cases for the entanglement of a single region. 

Let us now examine the property of I 3 for GHZ. Motivated by the above discussion of the 
EPR case, we note that in order to get a possibly nonzero I 3 in the GHZ case, we need some 
number of particles in all three regions. From the property that tracing out any number of 
particles leads to the same amount of entanglement, we conclude that 

S{A) = S{B) = S{C) = S{AB) = S{BC) = S{AC) . (B29) 
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From (B27) this implies 

h{A, B, C) = S{ABC) > 0. (B30) 

S{ABC) = 0, if there are no particles outside A, B, C, while S{ABC) ^ 0 (and equal to all 
the rest in (B29)), if there are particles in A, B, C and outside of these three regions. 

Note that GHZ states are very special, choosing a different multipartite entanglement 
pattern for the particles would generically lead to a non-definite sign of I 3 . 


c. Random pure state measure 


Consider a light cone which intersects with A and 5, and denote the area of its intersec¬ 
tions as uja,b respectively. From (4.16), assuming uja > ojbi the contribution from this light 
cone to I{A, B) is: 


I {A, B) D 2s (^6* ^B 

-\-6 ( cja + ~ ^ i 2 ^*^”^ ~ 


CJA + ~ 


(B31) 


From this expression one can check explicitly that the mutual information is non-extensive 
and monogamous, i.e. 


I 3 (A,B,C)< 0 , or I(A,B) + I(A,C) <I(A,BC). (B32) 

In contrast to previous two examples, I{A,B) is nonzero, if there exists a light cone, 
which intersects both A and B, and more than half of its area is inside AB. This implies 
that the mutual information stays zero for all times, if the separation of the regions is too 
large compared to their sizes. Also for large times we always get zero mutual information 
for compact regions. These results are reminiscent of holographic examples: holographic 
mutual information is monogamous [14], and the time dependence of mutual information 
shows similar qualitative behavior in holography, see e.g. [18]. 

In finite volume we can determine the saturation value of I(A,B) from (B23). For 
Va > Vb we get 

I(^A, 5) 2s 2 ^ystem^ 

fisystem^ ^ 

which has the same form as (B31). 


2 System -Vh 


Va + Vb ~ 2 System 


(B33) 


+0[Va + Vb-- 
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d. Example of two strips in d = 3 


As an example consider the mutual information for two parallel strips A and B of width 
a and length L, separated by a distance / in d = 3. In this case the RPS model will have 
zero mutual information unless j < ~ 0.207. If this holds the mutual information 

will be non zero starting at t = 1/ \/2, in contrast to the GHZ and EPR models which have 
non zero mutual information starting at t = //2 for all values of a//. In these latter models 
the mutual information decays asymptotically to zero with time, while in the RPS model is 
already zero for a finite time. This is illustrated in figure 30 where we plot the numerical 
evaluation of mutual information I {A, B) / (seqLa) as a function of time, r = t/a for these 
models. Note that the higher m GHZ models produce slower decay of mutual information, 
just as they give slower saturation for the entanglement entropy for one strip shown on 
Fig. 25. 


S(T)/(^eqLa) S(T)/(SeqLa) 



FIG. 30. Time evolution of mutual information divided Seq and per unit length for two strips of 
width a separated by a distance ^ in d = 3. The left panel is for a/Z = 1, while the right panel is 
for a/l — 0.2. The blue curve is the RPS model, and the GHZ curves for m — 1,2,3 are plotted in 
increasingly darker shades of orange. Note that for a/l = 0.2 (right panel) no mutual information 
is generated for the RPS model. 
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Appendix C: Explicit expressions for entanglement entropies after one scattering 


event 


( f) ( f) 

Here we give explicit expressions for S'jg and that we used in the discussion of the 
effect of a scattering event in Sec. VIA. 

For Si 3 we have 

'S'lS = -Pl logPl - (Pa - Pi) log(pa - pi) - P 2 logP 2 - (1 - Pa - P 2 ) log(l - Pa - P 2 ) , (Cl) 


where 


and 


Xi + a:2|Q;p 

(1 + |tt| 2)2 ’ 


Pa 


1 

1 + |tt |2 ’ 


P2 


2X3 + X4|q;| 

(1 + |tt| 2)2 ’ 


(C 2 ) 


X 2 = |(72ir + |(723l", *3 = |C/3l|"+|t/33|" , X, = |(741 r + |(743l" ■ (C3) 

Note that 

0 < a:i, 2 : 2 , 2 : 3 , 2:4 < 1 , + 2:2 + 2:3 + 2:4 = 2 . (C4) 

For I a I = 1 we explained in the main text that get a decrease in entropy. Indeed the above 
formulas simplify to 

- Pi^ log 

where we used that in this case pi < 1 / 2 . 

For we have 

SP = -pf logp/ - (1 - Pf) log(l - pf) (C 6 ) 

where 

— IcrP 

+ (i + |a|y <^‘ - = (TThV 

and 

0 < = |Cii|^ + iCisI^ ^1) 0 ^X 4 = < 1 . (C 8 ) 

Clearly for |q;| = 1, reduces to Sj^\ while for other values of a, one can always find 
values of xi and X 4 which either reduce or enhance it. 


c(/) _ 
>^13 ~ 


pilogpi + ( - 


Pl 


< 21 og 2 = S 


(i) 

13 ) 


(C5) 
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